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Abstract 

We present an overview of low-momentum two-nucleon and many-body interactions and their 
use in calculations of nuclei and infinite matter. The softening of phenomenological and effective 
field theory (EFT) potentials by renormalization group (RG) transformations that decouple low 
and high momenta leads to greatly enhanced convergence in few- and many-body systems while 
maintaining a decreasing hierarchy of many-body forces. This review surveys the RG-based tech- 
nology and results, discusses the connections to chiral EFT, and clarifies various misconceptions. 
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1 Introduction 

A new era is dawning for the theory of nuclear structure and reactions. Renewed interest in the physics 
of nuclei is fueled by experiments at rare isotope beam facilities, which open the door to new regions 
of exotic nuclei; by astrophysical observations and simulations of neutron stars and supernovae, which 
require controlled extrapolations of the equation of state of nucleonic matter in density, temperature, and 
proton fraction; and by studies of universal physics, which unite cold atom and dilute neutron physics [Ti- 
ll]. The interplay and coalescence of different threads: rapidly increasing computational power, effective 
field theory (EFT), and renormalization group (RG) transformations are enabling the development of 
new many-body methods and the revival of old ones to successfully attack these problems. 

A key to optimizing computational power for describing nuclei is a proper choice of degrees of 
freedom. While there is little doubt that quantum chromodynamics (QCD) is the correct underlying 
theory of strong interactions, the efficient low-energy degrees of freedom for nuclear structure are not 
quarks and gluons, but the colorless hadrons of traditional nuclear phenomenology. But this realization 
is not enough. For low-energy calculations to be computationally efficient (or even feasible in some 
cases) we need to exclude or, more generally, to decouple the high-energy degrees of freedom. 

Progress on the nuclear many-body problem has been slowed for decades because nucleon-nucleon 
(NN) potentials that reproduce elastic scattering phase shifts typically have strong short-range repulsion 
and strong short-range tensor forces. The consequence is substantial coupling to high-energy modes, 
which is manifested as highly correlated many-body wave functions and highly nonperturbative few- 
and many-body systems. In recent years, new approaches to nuclear forces grounded in RG ideas 
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Figure 1: NN phase shifts for the Argonne fig [H] (sohd), CD-Bonn fl9\ (dashed), and one of the chiral 
N^LO [20] (dotted) potentials in selected channels (using nonrelativistic kinematics). All agree with 
experiment up to about 300 MeV. 



and techniques have been developed. The RG allows continuous changes in "resolution" that decouple 
the troublesome high-momentum modes and can be used to evolve interactions to nuclear structure 
energy and momentum scales while preserving low-energy observables [SH?]. Such potentials, known 
generically as "low-momentum interactions," are more perturbative and generate much less correlated 
wave functions [8HT7]. This greatly simplifies the nuclear many-body problem, making structure and 
reaction calculations more convergent, while variations of the resolution provide new tools to assess 
theoretical errors. 

In this review, we survey the technical and phenomenological aspects of the low-momentum methods. 
Although there are multiple paths to low-momentum interactions, we focus on the RG-based techniques 
(known as "Viowfc" or "SRG" potentials), which provide new perspectives that mesh constructively with 
the developments of EFT for nuclear forces. When combined with advances in many-body methods and 
the increases in computer power, EFT and RG make feasible a controlled description, grounded in QCD 
symmetries, of nuclei across the nuclear many-body landscape. At the same time, the RG approach 
leads to reinterpretations of the physics or the role of different parts of the physics, such as what causes 
nuclear saturation. An unintended consequence is that many misconceptions or misinterpretations have 
arisen. A principal goal of this review is to address these. 



1.1 Nuclear forces 



Establishing an interparticle Hamiltonian, which is the most basic precursor to many-body calculations, 
is a difficult and on-going challenge for low-energy nuclear physics. The two-body sector has been 
"solved" in the sense that various interactions are available that reproduce phase shifts with x^/dof ~ 1 
in the elastic regime (up to roughly 300-350 MeV energy in the laboratory frame, see Fig. [1]). The 
unsettled frontier is three- and higher-body forces, although there remain important open questions 
about the systematic construction of NN potentials using EFT. 



Figure 2(a) shows nuclear interactions in the ^Sq channel for several phenomenological NN poten- 



tials. The longest range feature is one-pion exchange, which is justified by quantum chromodynamics 
(via the spontaneous breaking of chiral symmetry) and is a common feature of most potentials. The 
midrange part, which has a net attraction, has usually been associated with two-pion and/or heavy 
meson exchange (p, u, "cr"). The short-range part of the potentials in Fig. 2(a) is a repulsive core 
(often called a "hard core"). 

Nuclear structure calculations are complicated due to the coupling of low to high momenta by these 
potentials. This is made clear by the Fourier transform (that is, the Bessel transform in a given partial 
wave), as shown in Fig. 2(b) We feature the Argonne vis potential [18] because it is used in the most 
successful high precision (< 1% accuracy) nuclear structure calculations of nuclei with mass number 
A ^ 12 [22H2I]. For our purposes, the equivalent contour plot in Fig. 3(a) is a clearer representation 
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Figure 2: (a) Several phenomenological NN potentials in the ^Sq channel from Ref. [2Tj. (b) Momentum- 
space matrix elements of the Argonne vis (AV18) ^Sq potential after Fourier (Bessel) transformation.^ 



and we use such plots throughout this reviewQ The elastic regime for NN scattering corresponds to 
relative momenta k < 2fm~^. The strong low- to high- momentum coupling driven by the short-range 
repulsion is manifested in Fig. 3(a) by the large regions of non-zero off-diagonal matrix elements. A 
consequence is a suppression of probability in the relative wave function ("short-range correlations"), 



as seen for the deuteron in Fig. 3(b) 



The potentials in Fig. 2(a) are partial- wave local; that is, in each partial wave they are functions 
of the separation r alone. This condition, which simplifies certain types of numerical calculations^ 



constrains the radial dependence to be similar to Fig. 2(a) if the potential is to reproduce elastic phase 
shifts, and in particular necessitates a strong short-range repulsion in the S-waves. The similarity 
of all such potentials, perhaps combined with experience from the Coulomb potential, has led to the 
(often implicit) misconception that the nuclear potential must have this form. This prejudice has been 
reinforced recently by QCD lattice calculations that apparently validate a repulsive core |25H2B]- 

For finite-mass composite particles, locality is a feature we expect at long distances, but non-local 
interactions would be more natural at short distances. In fact, the potential at short range is far 
removed from an observable, and locality is imposed on potentials for convenience, not because of 
physical necessity. Recall that we are free to apply a short-range unitary transformation U to the 
Hamiltonian (and to other operators at the same time). 



1) 



and the physics described by H and H is indistinguishable by experiment. Thus there are an infinite 
number of equally valid potentials, and once we allow non-locality, a repulsive core and the strong low- 
to high-momentum coupling is no longer inevitable. 

The EFT approach uses this freedom to construct a systematic expansion of the Hamiltonian. A 
particular EFT is associated with a momentum scale Ab that is the dividing point between resolved. 



"'^In units where h ^ c = m = 1 (with nucleon mass to), the momentum-space potential is given in fm. In addition, we 
typically express momenta in fm^^ (the conversion to MeV is using he « 197MeVfm). 

^For example, in current implementations of Green's Function Monte Carlo (GFMC) calculations (32], the potential 
must be (almost) diagonal in coordinate space, such as the Argonne vis potential. 
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Figure 3: (a) Momentum-space matrix elements of the Argonne vis ^Sq potential, (b) Probability 
density for the S-wave part of the deuteron wave function for Argonne vis and smooth Viowfc potentials 
with several cutoffs A. The probability suppression at short distances is called "short-range correlation" . 



long-range physics, which is treated explicitly, and unresolved, short-range physics, which is expanded 
in contact interactions. Results are given order-by-order in Q/Ab, where Q is a generic momentum (or 
light mass) scale of the process being calculated. There is also a cutoff A needed to regulate the theory, 
which suppresses high momenta. Thus A acts as a resolution scale for the theory. If A is chosen to 
be less than Ab, then the truncation error for the EFT will be dominated by powers oi Q/A rather 
than Q/Ab. In principle one could take A as large as desired but in practice this only works if the 
renormalization and the numerics involved in matching to data are sufficiently under control [29j . 

In general, the forces between nucleons depend on the resolution scale A and are given by an 
effective theory for scale-dependent two-nucleon I^n(A) and corresponding many-nucleon interactions 
V3n(A), V4n(A) and so on [5|l30p31j . This scale dependence is analogous to the scale dependence of parton 
distribution functions. At very low momenta Q <^ m^r, the details of pion exchanges are not resolved 
and nuclear forces can be systematically expanded in contact interactions and their derivatives |5U] . 
The corresponding pionless EFT (for which Ab ~ m^) is very successful in capturing universal large 
scattering-length physics (with improvements by including effective range and higher-order terms) in 
dilute neutron matter and reactions at astrophysical energies [3011321435] . 

For most nuclei, the typical momenta are Q ~ and therefore pion exchanges are included 
explicitly in nuclear forces. The corresponding chiral EFT has been developed for over fifteen years as 
a systematic approach to nuclear interactions |30l[3Tl[36l[37]. This provides a unified approach to NN 
and many-body forces, and a pathway to direct connections with QCD through lattice calculations (see, 
for example, Ref. p8]). Examples of order- by-order improved calculations of observables are shown in 



Figs. 4(b) , 5(a) , and 5(b) However, some open questions remain [31]: understanding the power counting 
with singular pion exchanges [5UH1T] . including A degrees of freedom, the counting of relativistic 1/m 
corrections. Resolving these questions is important for improving the starting Hamiltonian for low- 
momentum interactions, but does not affect our discussion of RG technology. 

In chiral EFT [301 ED ESI EZ], the expansion in powers of Q/Ay, has roughly Ab < nip. As shown 



in Fig. 4(a) , at a given order this includes contributions from one- or multi-pion exchanges and from 
contact interactions, with scale-dependent short-range couplings that are fit to low-energy data for each 
A (experiment captures all short-range effects). There are natural sizes to many-body force contributions 
that are made manifest in the EFT power counting and which explain the phenomenological hierarchy 
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Figure 4: (a) Chiral EFT for nuclear forces, (b) Improvement in neutron-proton phase shifts shown 
by shaded bands from cutoff variation at NLO (dashed), N^LO (hght), and N^LO (dark) compared to 
extractions from experiment (points) [31]. The dashed hne is from the N^LO potential of Ref. |20] . 



of many-body forces. In addition, the EFT (extended to include chiral perturbation theory) provides 
a consistent theory for multi-pion and pion-nucleon systems and electroweak operators, as well as for 
hyperon-nucleon interactions [3111121^3] . 

The highest-order NN interactions available to date are at next-to-next-to-next-to-leading order, 
N^LO or (Q/Ab)^, for several different cutoffs (A = 450-600 MeV) and two different regulator schemes [20 



Representative results for NN phase shifts at NLO, N^LO, and N^LO are shown in Fig. 4(b) 
where error bands are determined by the spread in predictions for different A ~ Ab. Contour plots of 
momentum-space matrix elements for the softest, most commonly used N^LO potential and one with 
a higher cutoff are shown in Fig. |6l While they are much softer than Argonne f is in the ^Sq channel, 
there is still considerable off-diagonal strength above k = 2fm~^, which remainsproblematic for nuclear 
structure calculations (and the coupled ^Si— '^Di channel is generally worse) |j One might think the 
solution is to simply fit with a smaller A, but then the fit worsens significantly as the truncation error 
grows with Q/A. 

Instead we can use RG transformations to evolve to lower A while preserving the truncation error 
of the original Hamiltonian. We will roughly define "low-momentum interactions" as potentials that do 
not couple k <2 fm~^ to larger momenta. This is a relatively small (but significant) evolution for chiral 
potentials and a large one for phenomenological potentials. Low-momentum interactions are sometimes 
categorized as phenomenological interactions or regarded as an alternative to EFT interactions. Instead 



■^Note that the cutoff associated with the potential in Fig. 6(a) is A = 500 MeV, which might lead one to expect no 
strength above k « 2.5 fm^^. However, the regulator does not sharply cut off relative momenta. 
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Figure 5: (a) Differential cross section (in mb/sr) and vector analyzing power for elastic neutron- 
deuteron scattering at 10 MeV (top) and 65 MeV (bottom) at NLO (light) and N^LO (dark) from 
Ref. [36]. (b) Ground-state energy of ^Li at NLO and N^LO with bands corresponding to the A 
variation over 500-600 MeV compared to experiment (solid line, see Ref. [36] for details). 
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Figure 6: Two N^LO ^Sq potentials with (a) 500 MeV and (b) 600 MeV cutoffs [20]. 



they are an entire class of potentials associated with an initial Hamiltonian. There is no prejudice in 
the RG methods as to the starting Hamiltonian, but we favor the chiral EFT framework because of the 
consistent organization of many-body forces and operators and the capacity for systematic improve- 
ments. We emphasize that precision tests of all currently available Hamiltonians have only been made 
in few-body systems. Details of three-body forces and the four-body strength, for example, are quite 
uncertain o 

In the past, unitary transformations were used to soften NN potentials, which were then applied 
to calculations of few-body systems and nuclear matter. Results for observables depended on the 
transformations and this was often the context for discussing "off-shell" effects and which was the "true 
potential". From the modern perspective, this approach is misleading at best. These transformations 



''An exploratory calculation for ^He suggested a contribution of order 100 kcV to the binding energy from the long-range 
four-body force at N^LO, but with substantial uncertainty |45II46) . 
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Figure 7: (a) Nuclear many-body landscape and (b) degrees of freedom and corresponding scales in 
nuclei HHI. 



always lead to many-body interactions, even if absent in the initial Hamiltonianjf] While this fact 
was clearly recognized in past investigations [47], many-body forces were usually neglected, which is 
what led to the different results. In contrast, we account for three-body forces, but we also stress 
that maintaining exact equivalence is not required, because the original Hamiltonian already contains 
truncation errors. 



1.2 Nuclear structure challenges 



The nuclear many-body landscape is illustrated in Fig. 7(a) , where each square represents a nucleus. 
The challenge of contemporary nuclear structure physics is to describe this entire range of nuclei and 
beyond to neutron stars and supernovae in a controlled and unified way. The logarithmic scale in 
proton and neutron number emphasizes the different regions in A where different many-body methods 
are applicable. Also shown is the theoretical hierarchy that we propose to exploit: QCD — )■ chiral EFT 
low-momentum interactions, which then feed into (most of) the many-body methods, including some 
not listed. Within the broader challenge there are many more specific challenges, for example, mapping 
out the drip-lines (the boundaries of nuclear existence), the description of halo nuclei, the evolution of 
shell structure, the characteristics of pairing, and the description of nuclear decays and spontaneous 
fission. An intertwined challenge to controlled calculations of binding energies and spectra are nuclear 
reactions. This includes the interaction with external electroweak probes, which require the consistent 



^It is possible to find NN potentials that come close to the binding energy of "^H and '*He and other few-body nuclei 
without three-body forces. Beyond the problem of constructing consistent electroweak currents for these interactions, 
this entails a fine tuning to avoid the natural size of many-body forces, which may not persist when these interactions 
are applied elsewhere. 
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calculation of electroweak currents. 

These challenges are increased by the wide range of scales for strongly interacting particles, illus- 



trated schematically in Fig. 7(b) In particular, we highlight the weak binding of nuclei compared to 
typical QCD energy scales. For nuclear structure physics, we need methods appropriate to the energy 
and momentum scales of nuclei, which means isolating the relevant degrees of freedom. As we emphasize 
throughout this review, the decoupling of scales is at the heart of EFT and RG approaches, and leads 
us to the low-momentum interaction technology to be documented. 

Each of the many-body methods faces computational challenges, principally to the rapid growth 
of resources needed with increasing A. For methods such as the No-Core Shell Model and coupled- 
cluster methods, the key is to improve convergence within the spaces accessible. For heavier nuclei, 
density functional theory methods offer favorable scaling, but are challenged to reach reliable accuracy 
such as the sub-MeV level for binding energies. Improved accuracy is needed in the shell model and 
other methods to realize the use of nuclei as laboratories for fundamental symmetries, such as for 
isospin-symmetry-breaking corrections to superallowed decays, for neutrinoless double-beta decay, or 
for octupole enhancement factors of electric dipole moments. In all cases, we seek theoretical error 
estimates, particularly for extrapolations to systems where measurements will be limited or non-existent. 

Progress toward such controlled nuclear calculations has long been hindered by the difficulty of 
the nuclear many-body problem when conventional nuclear potentials are used. This has historically 
been accepted as an unavoidable reality. Indeed, conventional wisdom among nuclear physicists, as 
summarized by Bethe in his review of over 30 years ago [49j, holds that successful nuclear matter 
calculations must be highly nonperturbative in the potential. This is in contrast to the Coulomb many- 
body problem, for which Hartree-Fock is a useful starting point and (possibly resummed) many-body 
perturbation theory is an effective tool. The possibility of a soft potential providing a more perturbative 
solution to the nuclear matter problem was discarded at that time, and saturation firmly identified 
with the density dependence due to the tensor force [19]. Until recent RG-based calculations [9l[T6t 
[T7] . subsequent work on the nuclear matter problem [501453] had not significantly altered the general 
perspective or conclusions of Bethe's review (although the role of three-nucleon (3N) forces has been 
increasingly emphasized) . 

As already noted, nonperturbative behavior in the particle-particle channel for nuclear forces arises 
from several sources. First is the strong short-range repulsion, which requires at least a summation of 
particle-particle ladder diagrams jlD]. Second is the tensor force, for example, from pion exchanges, 
which is highly singular at short distances, and requires iteration in the triplet channels [5ll|55]. Third 
is the presence of low-energy bound states or nearly-bound states in the S-waves. These states imply 
poles in the scattering T matrix that render the perturbative Born series divergent. All of these 
nonperturbative features are present in conventional high-precision NN potentials. 

The philosophy behind the standard approach to nuclear matter is to attack these features head-on. 
This attitude was succinctly stated by Bethe [19]: 

"The theory must be such that it can deal with any NN force, including hard or 'soft' core, 
tensor forces, and other complications. It ought not to be necessary to tailor the NN force 
for the sake of making the computation of nuclear matter (or finite nuclei) easier, but the 
force should be chosen on the basis of NN experiments (and possibly subsidiary experimental 
evidence, like the binding energy of '^H)." 

In contrast, the EFT and RG perspective has a completely different underlying philosophy, which 
stresses that the potential is not an observable to be fixed from experiment (there is no "true potential"), 
but that an infinite number of potentials are capable of accurately describing low-energy physics |56j . 
In order to be predictive and systematic, an organization ("power counting") must be present to permit 
a truncation of possible terms in the potential. If a complete Hamiltonian is used (including many-body 
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Figure 8: Excitation energies (in MeV) in light nuclei calculated using the No-Core Shell Model (NCSM) 
with chiral EFT interactions (NN to N^LO and 3N to N^LO) compared to experiment [57]. Reprinted 
with permission from P. Navratil et al. copyright (2007) by the American Physical Society. 

forces), then all observables should be equivalent up to truncation errors. The EFT philosophy implies 
using this freedom to choose a convenient and efficient potential for the problems of interest. 

The use of energy-independent low-momentum interactions is a direct implementation of these ideas. 
Varying the cutoff can be used as a powerful tool to study the underlying physics scales, to evaluate 
the completeness of approximate calculations, and to estimate truncation errors from omitted higher- 
order contributions. These variable-cutoff interactions reveal the resolution or scale dependence of the 
first two sources of nonperturbative behavior, which are tamed as high momenta are decoupled. In 
free space, the third source of nonperturbative behavior remains independent of the cutoff because the 
pole positions of weakly and nearly bound states that necessitate fine tuning are physical observables. 
However, this fine tuning is eliminated in the medium at sufficiently high density. In short, a repulsive 
core is not constrained by phase shifts and is essentially removed by even a moderately low-momentum 



cutoff (note the A dependence in Fig. 3(b) ), the short-range tensor force is tamed by a sufficiently low 
cutoff, and the weakly and nearly bound states become perturbative as a result of Pauli blocking. For 
cutoffs around 2fm^^, which preserve phase shifts up to 330 MeV laboratory energy, the Born series in 
nuclear matter is well converged at second order in the potential, bringing the nuclear and Coulomb 
many-body problems closer together 

While evolving a soft potential from higher momentum is a new development in nuclear physics [5u58j, 
attempts to use soft potentials for nuclear matter were made in the mid sixties and early seventies [ITJ 
|59] . It had long been observed that a strongly repulsive core is not resolved until eight times nuclear 
saturation density [H]. Thus, saturation is not driven by a hard core (unlike liquid ^He). However, these 
soft potentials were abandoned because they seemed incapable of quantitatively reproducing nuclear 
matter properties. Their requiem was given by Bethe [49J: 

"Very soft potentials must be excluded because they do not give saturation; they give too 
much binding and too high density. In particular, a substantial tensor force is required." 

From the EFT perspective, a failure to reproduce nuclear matter observables should not be interpreted 
as showing that the low-energy potential is wrong, but that it is incomplete. This misconception still 
persists and has led to the conclusion that low-momentum NN interactions are "wrong" because they 
do not give saturation in nuclear matter and finite nuclei are overbound for lower cutoffs. The missing 
physics that invalidates this conclusion is many-body forces. 

In a low-energy effective theory, many-body forces are inevitable; the relevant question is how large 
they are. It is established beyond doubt that 3N forces are required to describe light nuclei [22l - |2H [57 t 
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Figure 9: Schematic illustration of two types of RG evolution for NN potentials in momentum space: 
(a) Viowk running in A, and (b) SRG running in A. At each Aj or Aj, the matrix elements outside of the 
corresponding lines are zero, so that high- and low-momentum states are decoupled. 

[60l|6T], as shown, for example, in Fig. |8l For variable-cutoff potentials, three-body (and higher-body) 
interactions evolve naturally with the resolution scale. 

1.3 Renormalization group approaches 

A fundamental tenet of renormalization theory is that the relevant details of high-energy physics for 
calculating low-energy observables can be captured in the scale-dependent coefficients of operators 
in a low-energy Hamiltonian [29]. This principle does not mean that high-energy and low-energy 
physics is automatically decoupled in every effective theory. In fact, it implies that we can include as 
much irrelevant coupling to incorrect high-energy physics as we want by using a large cutoff, with no 
consequence to low-energy predictions (assuming we can calculate accurately). But this freedom also 
offers the possibility of decoupling, which makes practical calculations more tractable by restricting 
the necessary degrees of freedom. This decoupling can be efficiently achieved by evolving nuclear 
interactions using RG transformations designed to handle similar problems in relativistic field theories 
and critical phenomena in condensed matter systems]^ 

The general purpose of the RG when dealing with the large range of scales in physical systems was 
eloquently explained by David Gross p5] : 

"At each scale, we have different degrees of freedom and different dynamics. Physics at a 
larger scale (largely) decouples from the physics at a smaller scale. . . . Thus, a theory at a 
larger scale remembers only finitely many parameters from the theories at smaller scales, 
and throws the rest of the details away. More precisely, when we pass from a smaller scale 
to a larger scale, we average over irrelevant degrees of freedom. . . . The general aim of the 
RG method is to explain how this decoupling takes place and why exactly information is 
transmitted from scale to scale through finitely many parameters." 

The common features of RG for critical phenomena and high-energy scattering are discussed by Steven 
Weinberg in an essay in Ref. ^64j. He summarizes: 




"The method in its most general form can I think be understood as a way to arrange in 
various theories that the degrees of freedom that you're talking about are the relevant degrees 
of freedom for the problem at hand." 

^For an early discussion of decoupling based on Okubo unitary transformations, see Ref. [52] . 
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Figure 10: Two types of RG evolution applied to one of the chiral N^LO NN potentials (550/600 MeV) 
of Ref. [H] in the ^Si channel: (a) V^owfc running in A, and (b) SRG running in A (see Fig. [27] for plots 
in k"^, which show the diagonal width of order A^). 



This is the heart of what is done with low-momentum interaction approaches: arrange for the degrees 
of freedom for nuclear structure to be the relevant ones. This does not mean that other degrees of 
freedom cannot be used, but to again quote Weinberg [6l]: "You can use any degrees of freedom you 
want, but if you use the wrong ones, you'll be sorry." 

There are two major classes of RG transformations used to construct low-momentum interactions, 
which are illustrated schematically in Fig. [91 In the V^owfc approach, decoupling is achieved by lowering 
a momentum cutoff A above which matrix elements go to zero. In the SRG approach, decoupling is 
achieved by lowering a cutoff A (in energy differences A^) using flow equations, which means evolving 
toward the diagonal in momentum space. The technology for carrying these out is outlined in Section [3l 



but the effects can be readily seen in the series of contour plots in Figs. 10(a) and 10(b) 



With either approach, lowering the cutoff leaves low-energy observables unchanged by construction, 
but shifts contributions between the interaction strengths and the sums over intermediate states in loop 
integrals. The evolution of phenomenological or chiral EFT interactions to lower resolution is beneficial 
because these shifts can weaken or largely eliminate sources of nonperturbative behavior, and because 
lower cutoffs require smaller bases in many-body calculations, leading to improved convergence for 
nuclei. The RG cutoff variation estimates theoretical uncertainties due to higher-order contributions, 
to neglected many-body interactions or to an incomplete many-body treatment. When initialized with 
different orders of chiral EFT interactions, we have a powerful tool for extrapolations to the extremes 
and for assessing the uncertainties of key matrix elements needed in fundamental symmetry tests. 

The idea of effective interactions in a limited model space is an old and well-exploited one in nuclear 
physics. However, we will emphasize the flexibility of the RG compared to effective interaction methods. 
The continuous "cutoff" variation (in quotes because it may not be an explicit cutoff) is a valuable new 
tool for nuclear physics (see Section I2l6l) . The RG methods are versatile and suggest new ways to make 
progress (for example, using the in-medium SRG, discussed in Section 1^^ . In addition, RG combined 
with EFT is a natural framework for uncovering universal behavior. 

We note that the RG is an integral part of any EFT. Matching of the EFT at a given truncation 
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level (to data or to an underlying theory) but at different regulator cutoffs establishes the RG evolution 
(or "running") of the EFT couplings. This includes the shift of strength between loop integrals and 
couplings and between two and many-body interactions. However, because the EFT basis is truncated, 
the error at the initial cutoff is not preserved with the running, in contrast to the momentum-space RG 
evolution using V^owfc or SRG techniques, which keep all orders. 

1.4 Scope of the review 

This review focuses on the use of RG methods to derive low-momentum interactions. There are other 
approaches to softened interactions that share many of the same advantages and issues, which we do not 
address in detail. Fortunately there are recent reviews to cover these omissions, for example, Ref. |3T] 
on nuclear forces, Refs. [65] and [66] on ab initio and shell-model applications, and Ref. [67] on density 
functional theory. 

In Section [2l we discuss general principles of the RG as applied to nuclear forces. The various tech- 
nologies for evolving NN potentials are outlined in Section [3l Many-body interactions and operators are 
addressed in Section |H Sections [5] and [6] review many-body advances with applications to infinite matter 
and finite nuclei. We conclude in Section [7] with a summary of the main points, through a consideration 
of misconceptions and clarifications, on-going developments, and important open questions. 



2 Renormalization group: motivation and principles 
2.1 Decoupling 

High-momentum degrees of freedom do not automatically decouple from low-energy observables, espe- 
cially for NN potentials that have significant high-momentum off-diagonal strength (typically beyond 
k > 2fm~^). The coupling between low- and high-momentum states is a consequence of quantum 
fluctuations, as manifested in intermediate state summations in perturbation theory, e.g., 

hi £a 
q=0 ^ 

This coupling is regulated by the cutoff or resolution scale A that is present in all NN interactions. 

A central theme of this review is that nonperturbative features of nuclear forces can be radically 
altered by explicitly decoupling low- and high-momentum degrees of freedom. As we will see, many-body 
calculations become more perturbative, basis expansions converge more rapidly, and strong short-range 
correlations are smeared out, rendering variational methods more effective. However, if NN potentials 
at a large scale A are simply truncated at some fcmax <^ A, the calculated low-energy observables 
(such as phase shifts and the deuteron binding energy) are drastically changed, as shown in Figs. [H] 
and [121 Such a brute-force cutoff does not disentangle high-energy (short-range) features from low- 
energy (long-range) observables. To do so, it is necessary to integrate out (and thus decouple) irrelevant 
high-momentum details from their effects on low-energy observables. This is achieved by the two classes 
of RG transformations illustrated in Fig. [9] and described in greater detail in Section [3l 

A practical test for decoupling is whether changing high-momentum matrix elements of a given 
potential changes low-energy observables. The SRG allows for particularly convincing demonstrations of 
decoupling, because the unitary evolution of the Hamiltonian preserves the deuteron binding energy and 
phase shifts for all energies. The SRG-evolved potential Vg explicitly decouples high-energy dynamics 
from low-energy observables, which means that we can set to zero the high-momentum parts (so that 
we have a potential like V^owfc) with no significant changes in the low-energy physics. 
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Figure 11: NN pliase sliifts for tlie Argonne fis potential [18] and the SRG-evolved potential for 
A = 2fm~^ with all momenta included (two indistinguishable solid lines) and with the exclusion of 
momenta k > k^s.x = 2.2 fm^^ (Argonne fig dotted, SRG-evolved dashed) [68] . 
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Figure 12: Expectation values in the deuteron of the kinetic (T) (dashed), potential {V) (dot-dashed), 
and total energy (solid) evaluated in momentum space as a function of the maximum momentum fcmax, 
see Eq. ([3]). Results are shown for the Argonne fig potential [18] (left), the SRG-evolved potential Vg 
for A = 2fm~^ (middle), and the smooth-cutoff Vio^k interaction with A = 2fm~^ (right) [68] . 



In Ref. [68], the decoupling of high-energy details from low-energy phase shifts and the deuteron 
binding energy was demonstrated by setting Vs{k, k') to zero for all k, k' above a specified momentum 
kma.x (using a smooth regulator function). Phase shifts for fcmax = 2.2 fm~^ are shown in Fig. [TT] for 
the initial Argonne Vis potential and for the SRG-evolved Vs with A = 2fm~^. The phase shifts for 
the initial potential in the lower partial waves bear no relation to the result without a fcmax cutoff. In 
contrast, the low-energy phase shifts for the SRG-evolved potential are unchanged, even though the 
high-energy phase shifts above /cmax are now zero. 

The deuteron binding energy provides another clear example of how the contributions of different 
momentum components to a low-energy observable depend on the resolution scale (as measured by A 
or A, see Fig. In Fig. [T^l we show the kinetic, potential, and total energy from an integration in 
momentum space including momenta up to fcmax- That is, we plot 

^d(A; < A;^ax) = / dk dk>J(k;A)(A;253(k-k') + K(k,k'))^d(k';A), (3) 

Jo Jo 

where ipdik; A) is the momentum-space deuteron wave function from the corresponding potential Vg 
(without /cmax)- Figure [12] shows that if one excludes momenta greater than 2fm~^ in the Argonne vis 
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Figure 13: Ground-state energy of as a function of the maximum momentum fcmax for three different 
values of A [H]. The cutoff function is exp[— (fc^//;;^^^)"] with n = 8. The initial potentials are (a) the 
N^LO NN potential of Ref. [20] and (b) the Argonne Vis potential [18]. 



deuteron wave function, the deuteron is 9.9 MeV unbound (that is, the integrated kinetic energy up to 
2fm~^ is 11.5 MeV while the potential energy is —1.6 MeV). In contrast, using Vs with A = 2fm~^, one 
sees that the converged result is dominated by contributions from much lower momenta. Note that the 
Vs potential has no appreciable contributions above A, even though the near-diagonal matrix elements 
of the potential Vs{k, k') for /c, k' > fcmax are sizable. This again validates decoupling. It is also evident 
that the Vs and Viowfc results are very similar for A ~ A, where A is the momentum cutoff for V^owfc- 
Lastly, we mention that the decoupling carries over to few-body systems [HES], as shown similarly for 
the ground-state energy in Fig. [131 Note that the different converged energies for each A reflects the 
truncation of the SRG at the NN level, see Section 12.51 



2.2 Connections to effective field theory 

The RG plays a central role in EFT. It determines the running of low-energy couplings (see, for example, 
Ref. [30]) and the RG scaling of operators can be used to identify power counting schemes [70l[7l]. In 
this section, we discuss how the nonperturbative RG evolution works in the context of chiral EFT inter- 
actions. In EFT, it is optimal to fit the low-energy couplings at larger cutoffs, because this minimizes 
the truncation errors that can scale as rapidly as (Q/A)", and because larger cutoffs include maximal 
long-distance physics up to the given order in the EFT expansion. This is illustrated in Fig. [T3] by 
comparing an N^LO potential with a larger cutoff A = 500 MeV (2.5 fm"^) to the sharp cutoff N^LOW 
potential of Ref. [72] with A = 2.0fm"\ At the same EFT order, the fit for a lower cutoff requires 
more fine-tuning and leads to a potential that is not very smooth in momentum space. This indicates 
a breakdown of the gradient expansion. 

The RG evolution of EFT interactions to lower cutoffs is beneficial because the resulting low- 
momentum interactions preserve the truncation error of the initial potential and therefore include 
the advantages of fitting at larger cutoffs, while they have lower cutoffs that are advantageous for 
applications to nuclear structure (see Sections |5] and |6]). As shown in Fig. [TU the evolved interactions 
are smooth without the strong momentum dependencies observed for the sharp cutoff N^LOW potential 
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Figure 14: Momentum-space matrix elements in the ^Si channel for (a) the N'^LO potential (A = 
500 Me V) of Ref. [20], (b) the sharp cutoff N^LOW potential (A = 2.0 fm'^) of Ref. [72], and (c) for 
the N^LO potential (A = 500 Me V) evolved by a smooth Viowfc to A = 2.0fm"\ 
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Figure 15: Flow of Viowfc(^' = 0, /c = 0;A) compared to the corresponding momentum-independent 
contact interaction Co(A) at LO and NLO, where this coupling is determined entirely from RG invariance 
and fits to the scattering length (at LO) plus effective range (at NLO) [73] . 



at the same cutoff scale, and the NN phase shifts are reproduced with the same accuracy. 

We can gain further insights into the interplay of the RG and EFT by considering chiral EFT as 
providing a general operator basis that can be used to expand the RG evolution. At a given order 
(Q/Ab)", chiral EFT includes contributions from one- or multi-pion exchanges and from contact inter- 
actions, with short-range couplings that depend on the resolution or cutoff scale. As part of the RG 
evolution, short-range couplings included in the initial potential evolve. This is illustrated in Fig. [15] by 
comparing the flow of \4owfc(^' = 0, A; = 0; A) with the corresponding momentum-independent contact 
interaction Co(A) in subsequent orders of pionless EFT. In addition, the RG generates higher-order 
short-range contact interactions so that observables are exactly reproduced and the truncation error is 
unchanged. Consequently, the cutoff variation can be used to estimate theoretical uncertainties due to 
higher-order short-range many-body interactions. We will discuss using cutoff dependence as a tool in 
Section [2S1 

For cutoffs large compared to the pion mass, the long-range parts in the initial chiral potentials 
are preserved by the RG evolution. This is demonstrated in Fig. 16(a) for deuteron wave functions at 
different resolution scales. While the short-range behavior changes with the cutoff, the long-range parts 
are governed by low-energy pion physics and remain unchanged. The results discussed in Section H] for 
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Figure 16: (a) Deuteron wave functions in coordinate space for smooth Viowk interactions at several 
different cutoffs using an exponential regulator with nexp = 8. The initial interaction is the N^LO 
potential of Ref . [20] • (b) Matrix element of the bare quadrupole moment operator as a function of the 
cutoff for several smooth V^owfc regulators. For details see Ref. [6]. 



many-body interactions and operators can be explained based on the above observations. 



2.3 Universality 

Low-momentum interactions starting from different initial potentials are found to be quantitatively 
similar [5HZ]. We use the terminology universality in this context (independent of the high-momentum 
details in the initial potential) and not in the sense of critical phenomena (independent of the system 
details). Figures [T7] and (TS] show various phenomenological NN potentials and different N^LO NN 
interactions of Entem and Machleidt (EM) [20] and of Epelbaum et al. (EGM) ^3]. Each potential 



accurately reproduces low-energy NN scattering, as shown in Fig. 4(b) When the potentials are evolved 
to lower cutoffs, the resulting low-momentum interactions become universal. In addition, the RG 
evolution weakens the off-diagonal coupling between low and high momenta. Section [3] reviews different 
low-momentum technologies for achieving this decoupling. 

A similar universal behavior and decoupling is found for the low-momentum parts of SRG evolved 



interactions, see Fig. 28(a) We attribute the collapse to universal low- momentum interactions to the 
common long-range pion physics in the initial potentials, and to a similar description of low-energy NN 
observables up to the resolution scale, while the high-momentum parts are decoupled. This is similar to 
the inverse scattering problem, where the theory is restricted to the energy range set by the resolution 
scale. Important open problems are to understand the observed universality in terms of an RG scaling 
analysis, to identify eigenoperators of the RG evolution, and to extend the study of universality to 3N 
forces. 



2.4 Perturbativeness 

Nuclear many-body calculations are complicated by strong short-range repulsion and strong short-range 
tensor forces found in most NN potential models. However, both of these sources of non-perturbative 
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Figure 17: Diagonal (left) and off-diagonal (right) momentum-space matrix elements for various phe- 
nomenological NN potentials initially (upper figures) and after RG evolution to low-momentum inter- 
actions Vxawk ^i6j (lower figures) for a smooth regulator with A = 2.0 fm~^ and n, 
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Figure 18: Diagonal (left) and off-diagonal (right) momentum-space matrix elements of different N'^LO 
NN interactions (EM [20] and EGM [H]) initially (upper figures) and after RG evolution to low- 

lowfc [HIE] (lower figures) for a smooth regulator with A = 2.0 fm~^ and ngxp = 4. 
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Figure 19: (a) Particle-particle contributions to the energy per nucleon in symmetric nuclear matter 
as a function of the Fermi momentum fcp for the initial Argonne fig potential and the RG-evolved 
Viowk with A = 2.1 fm~ [9]. (b) Pair-distribution function g(r) in nuclear matter for kp = 1.35 fm~ at 
different resolutions, for details see Ref. [lOj . 



behavior are resolution-dependent because they depend on the degree of coupling between low- and 
high-momentum states [6l[71[9l[T2]. Consequently, KG methods can be used to improve perturbative 
convergence and reduce the short-range strength of the associated correlations in the wave functions. 



as shown for symmetric nuclear matter in Figs. 19(a) and 19(b), respectively. 



VK) = r],iE)\T,) , (5) 



We can quantify the perturbativeness of the potential as we evolve to lower A for Viowfc (or A for 
SRG) interactions by using the eigenvalue analysis introduced long ago by Weinberg [75J and applied 
to Viowfc and SRG potentials in Refs. l6l[7l|9l[T2] . Consider the Born series for the T matrix at energy E 
with Hamiltonian H = Hq + V , 

TiE) = V + V^^V + .... (4) 

By finding the eigenvalues and eigenvectors of the operator {E — Hq)^^V, 

1 

and then acting with T{E) on the eigenvectors, 

r(E)|r,) = (1 + vAE) + vliE) + . . .) v|r.) , (6) 

it follows that nonperturbative behavior at energy E is signaled by one or more eigenvalues with 
\ri^{E)\ ^ 1 [75j. A rearrangement of Eq. ([5]) gives a simple interpretation of the eigenvalue rj^{E) as 
an energy- dependent coupling that must divide V to produce a solution to the Schrodinger equation at 
energy E. For negative energies, a purely attractive V gives positive real ri^{E) values, while a purely 
repulsive V gives negative eigenvalues. For this reason, we refer to negative eigenvalues as repulsive 
and positive ones as attractive, although the eigenvalues become complex for positive E. 



^This is in contrast to non-perturbative features like low-energy bound or nearly bound states in the S-waves and the 
pairing instability at finite density that are insensitive to the short-distance details. 
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Figure 20: Trajectories of the largest repulsive Weinberg eigenvalues in the (a) ^5*0 and (b) ^Si-'^Di 
channels as a function of energy for V^iowfc evolved from the Argonne fig potential [12]. The results for 
selected cutoffs are indicated by the different symbols. The positions of the symbols on each trajectory 
mark the eigenvalues for center-of-mass energies Ecm = 0,25,66,100 and 150 MeV, starting from the 
filled symbol at OMeV. The trajectory with stars are eigenvalues for the N^LO potential of Ref. [20] . 




A (fm"') A (fm"') 

(a) (b) 

Figure 21: Largest repulsive Weinberg eigenvalues for ii^ = in the (a) ^Sq and (b) '^Si-'^Di channels 
as a function of cutoff for Viowfc evolved from chiral EFT interactions 0. Results are shown for the 
N^LO potential of Entem and Machleidt [20], for N^LO potentials of Epelbaum et al. [44] with different 
cutoffs A/A (as indicated in MeV), and for an N^LO potential [71]. For comparison, we have plotted 
the largest repulsive Weinberg eigenvalues for Viowfc evolved from the Argonne vis potential. 



Figure [201 shows the trajectories of the largest repulsive Weinberg eigenvalue in the ^Sq and ^Si- 
^Di channels as a function of (positive) energy for Viowfc interactions with various cutoffs evolved from 
the Argonne fig potential. The magnitude of the largest repulsive eigenvalue at all energies decreases 
rapidly as the cutoff is lowered. This reflects the decrease of the short-range repulsion present in the 
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Figure 22: (a) Evolution of the two largest Weinberg eigenvalues with A in the ^Si-^Di channel in free- 
space (solid) and at saturation density (dashed). The lighter curves correspond to attractive eigenvalues 
and the darker curves to repulsive eigenvalues, (b) Dependence on density of the two largest Weinberg 
eigenvalues in the ^Si-'^Di channel. Results are based on Viowfc evolved from the Argonne fig potential. 
For details see Ref. M- 



initial potential. In the ^Sq channel, the trajectory lies completely inside the shaded unit circle for 
cutoffs near 4fm~ and below, which implies that the Born series becomes perturbative with respect 
to the repulsive part of the potential. In the '^Si-^Di channel, the largest repulsive eigenvalues are just 
within the unit circle for A = 3fm~^. The decrease from A = 3fm~^ to 2fm~^ is significant for ensuring 
a convergence of particle-particle ladders in nuclear matter, because second-order short-range tensor 
contributions peak at intermediate-state momenta k !^ 2.5 — 3.5 fm~^ in nuclear matter [76] . 

Chiral EFT interactions typically have cutoffs below 3.5 fm~^ and therefore are expected to be soft 
potentials. However, as shown in Fig. the N^LO potentials of Epelbaum et al. [44] have substantial 
repulsive Weinberg eigenvalues in both S-waves. These unexpected features were traced to the singular 
central and tensor interactions (at N^LO) in Ref. [12]. The results in Fig. |2T] show there is a major 
decrease in the largest repulsive eigenvalues as the cutoff is run down from A = 3fm~^ to 2fm~^, which 
shows that it is advantageous to evolve chiral EFT interactions to lower cutoffs using the RG. 

The systematics of the Weinberg eigenvalues in free space and at finite density are compared for 



the ^Si-^Di channel in Fig. 22(a), which shows the evolution with A of the two largest eigenvalues 
evaluated at the deuteron pole, E = Bd- In free space, we observe that the repulsive eigenvalue is 
initially large but decreases rapidly with cutoff, while the large attractive eigenvalue corresponding to 
the deuteron pole remains invariant, ri^{Bd) = 1. The behavior is similar at finite density, but now the 
attractive eigenvalue is tamed by Pauli-blocking effects. This result is general at sufficient density, as 
shown in Fig. 22(b) We therefore reach a promising conclusion: Large-cutoff sources of nonperturbative 
behavior can be eliminated using the RG, while physical sources due to weakly and nearly bound states 
are suppressed in the medium by Pauli-blocking. These results are qualitatively unchanged with the 
inclusion of low-momentum 3N interactions [9l[16] (see Section |5]). Moreover, the Weinberg analysis has 
been applied to pairing in Ref. [77] . 

In addition to improving perturbative convergence, RG-evolved potentials result in more effective 
variational calculations [11] and rapidly converging basis expansions [lll[T5] , as discussed in Section |6l 
Finally, we mention that as a by-product of the Weinberg analysis in Ref. [12], it was found that 
RG-evolved interactions can be accurately described by low-rank separable expansions, a result that 
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Figure 23: Eliminating degrees of freedom leads to three-body forces. 



allows significant simplifications in applications ranging from few-body scattering [78j to the treatment 
of pairing in finite nuclei [701120] • 

2.5 Many-body interactions and operators 

In a low-energy effective theory of finite-mass composite particles such as atoms or nucleons, three-body 
(or higher-body) interactions are defined as a contribution to the Hamiltonian that is not accounted 
for by the sum of pairwise interactions. The polarization of interacting atoms or molecules provides an 
intuitive example. The Axilrod- Teller potential is a three-body version of the familiar two-body long- 
range van der Waals force between atoms; its physical origin is triple-dipole mutual polarization [81] . 
Because it contributes at third order in perturbation theory and the fine structure constant is small, there 
is a rapidly decreasing hierarchy of many-body forces. Indeed, this contribution is usually negligible 
in metals and semiconductors although not in rare gas solids [S2]. For solid xenon, its contribution 
is calculated to be about 10% of the ground-state energy [S3], comparable to the typical 3N force 
contribution to the binding energy. Note that polarization with additional atoms means that the 
Hamiltonian for an A-hodj system will inevitably lead to A-hodj forces. 

The EFT perspective (see Section 11.11) confirms that operators (including the Hamiltonian) in any 
low-energy effective theory, if it is systematic rather than a model, will have many-body components. 
The generic origin of these many-body operators is a restriction of degrees of freedom, such as 3N forces 
arising from the elimination of the A or anti-nucleon components, as illustrated in Fig. [231 Integrating 
out or decoupling high momentum modes with the RG is just another example. The EFT expansion 
systematically includes all such contributions as a combination of long-range (for example, due to pion 
exchanges) and short-range (contact) terms as on the right side of Fig. [231 even when the origin of 
the short-range parts is unclear. Power counting in chiral EFT establishes a decreasing hierarchy of 



many-body operators (see Fig. 4(a) for many-body forces) that permits truncation at a tractable level 
(which in present nuclear structure calculations means 3N forces and two-body current operators). 

The strength of many-body components in an initial operator will shift with any change in how high- 
energy degrees of freedom are coupled to the low-energy degrees of freedom; particular examples are the 
running of a V^owfc cutoff A or an SRG flow parameter A. The correlation plot in Fig. 24(a) of ^H and ^He 



binding energies calculated using Viowfc NN-only interactions shows the change in the (omitted) many- 
body contribution with A. This reproduces the empirical Tjon line from phenomenological potentials. 
The technical challenge is to carry out the evolution of many-body forces and operators in an RG 
implementation. The physics challenge is to establish that the EFT hierarchy of many-body components 
is not affected by the evolution, so that a tractable truncation is still possible. This typically means 
that the evolution is not extended below A or A of about 1.5 fm~ . (For cutoffs below the pion mass, 
three-body forces will increase to leading order in pionless EFT [301 [32].) 

The consistent RG evolution of Viowk many-body interactions has not yet been achieved. The 
underlying difficulty is that the technology used to construct requires the solution of the full 

3N problem (bound state wave functions plus scattering wave functions in all breakup channels) to 
consistently evolve 3N forces (without simplifications). In V^owA; calculations to date, the three-body 
evolution is therefore approximated by fitting the leading chiral EFT 3N forces at each cutoff while 
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(a) (b) 

Figure 24: Correlation of and ^He binding energies, (a) The cutoff dependence obtained from low- 
momentum NN interactions Viowk compared to the Tjon line [8]. (b) This correlation is driven by large 
scattering lengths, as demonstrated by the band obtained in pionless EFT [HI]. 



evolving the two-body interaction exactly [SKin]. This takes advantage of the EFT expansion being a 
complete operator basis for 3N forces, so that the leading effects of the evolution are simply a change 
in the operator coefficients (see Section W7L\ for more details). 

The SRG offers a new path to the consistent running of many-body interactions and operators 
because they evolve through unitary transformations that require only a representation in a convenient 
basis. This has recently been demonstrated in practice |85]; we discuss the details in Sections 14.21 and 
14. 4[ To see how the two-, three-, and higher-body potentials are identified and evolved, it is useful to 
decompose the running SRG HamiltoniaiJ^ Hx in second-quantized form. Schematically (suppressing 
indices and sums), 

Hx = (T) a^a + (l/f ^) a^a^aa + (vl^^) a^a^a^aaa + . . . , (7) 

where a\ a are creation and annihilation operators with respect to the vacuum in a single-particle basis. 
This defines (T), (V^^^''), (vlf^), ...as the one-body (kinetic energy), two-body, three-body, ...matrix 
elements at each A. The SRG evolution is dictated by commutators involving Hx (see Eq. ([21])); when 
they are evaluated using Hx from Eq. ([7]), we see that even if initially there are only two-body potentials, 
higher-body interactions are generated with each step in A. Thus, when applied in an A-body subspace, 
the SRG will "induce" A-hodj forces. But we also see that (T) is fixed, (V";^ ) is determined only in the 
A = 2 subspace with no dependence on {v!f^), {v!f^) is determined in A = 3 by subtraction for given 

(2) 

{Vx ), and so on. With this formulation, Hx is a free-space Hamiltonian, independent of A. It is also 
possible to normal order in the medium, which changes the definition of the creation and annihilation 
in Eq. ([7]) and shifts higher-body pieces to the zero-, one-, and two-body levels (see Section IT3]) . 

Although one often characterizes the changes with A or A by saying that the RG evolution "induces" 
many-body forces, this should not be interpreted as meaning the initial Hamiltonian has a "true" three- 
body force. The initial and evolved interactions are equally effective interactions and if the EFT 
truncation error is preserved in the evolution, there is no physics reason to prefer one over the other. 
(Indeed, referring to the initial interaction as the "bare" interaction is misleading in the same way.) 
Each Hamiltonian has different associated many-body components and operators, although the long- 
range parts are the same, as they are not modified by the evolution. A consequence of the latter is that 
operators that probe only long distances will not evolve, as seen by the weak cutoff dependence until 

^The SRG Hamiltonian is denoted equivalently by Hx or Hg, where A = l/s^/"*. 
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small A for matrix elements of the bare quadrupole moment operator in Fig. 



16(b) 



2.6 Cutoff dependence as a tool 

The RG equations for Viowfc or SRG interactions start with a Hamiltonian as an initial condition, which 
together with the corresponding operators fixes the values of observables. In principle there should be 
no change in observables as the Hamiltonian is evolved. (More precisely, in the usual Viowfc formulation, 
low-energy observables are unchanged while for the SRG all observables are preserved.) In practice 
there will be approximations both in the implementation of the RG and then in the subsequent calcu- 
lations of nuclear structure observables. That is, cutoff dependence arises because of the truncation (or 
approximation) of "induced" many-body forces (see Section 12. 5p or because of many-body approxima- 
tions. Because observables should be by construction independent of the RG cutoff or flow parameter, 
we can use changes as a diagnostic of approximations and to estimate theoretical errors. 

In more general applications of the RG, it is expected that significant changes in the cutoff are needed 
to glean useful information. In the examples considered here, however, the range of cutoff variation 
is a factor of two or in some cases even smaller. However, even small changes in the RG parameter 
can make dramatic shifts in the physics of nuclear structure observables, such as integrating out the 
short-range repulsion or short-range tensor forces. 

Here are some examples of using RG cutoff dependence as a tool: 



Cutoff dependence at different orders in an EFT expansion. In Figs. 4(b), 5(a), and 5(b), the 
spread of predictions for EFT interactions with a range of cutoffs indicates whether the results 
improve as expected. This carries over directly to low-momentum interactions constructed from 
initial EFT Hamiltonians at different orders. 

Running of ground-state energies with cutoff in few-body systems. The variation in Fig. 135) shows 
the scale of the (net) omitted many-body forces, which is then an estimate of the error from this 
approximation. This allows a comparison with expectations from EFT power counting. 

Tjon line and the nature of correlations (see Fig. |2^ . The cutoff dependence obtained from 



low-momentum NN interactions Viowfc in Fig. 24(a) explains the empirical (solid) Tjon line [S]. 



This correlation is driven by large scattering lengths, as demonstrated by the band in Fig. 24(b) 
obtained in pionless EFT [8^ . 

• Calculations of nuclear matter. The decrease in cutoff dependence with improved approximations 
(see Fig. HOj) helps validate the many-body convergence and establishes a lower bound on the error 
from short-range many-body interactions. 

• Calculations of nuclei. For example, the cutoff variation in helium halo nuclei with NN-only 
interactions identifies sensitivity to 3N force effects [I5l|86] (see Fig. HOj). We refer to Section [Q] 
for more examples of cutoff dependence used as a diagnostic of missing many-body forces. 

• Identification of (non-)observables. If a quantity changes more rapidly with the cutoff than any 
possible error, it is not an observable. An example is the D-state probability of the deuteron, 
whose cutoff dependence as shown in Fig. \57\ demonstrates it is not an observable. In contrast, 
the independence of the asymptotic D/S-state ratio is evidence it is an observable. 

• Cutoff dependence can be used in EFT as a tool to determine power counting (see Ref. [7T] and 
references therein). This entails running a cutoff A lower, rescaling in units of A, and looking for 
fixed points of the RG flow; that is, for flows that end in scale independence. The EFT is then 
expanded around such a fixed point. Examples for the nuclear case are described in Ref. [71] . 

The common lesson here is that it is advantageous to do calculations for a range of cutoff values. 
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3 Low- momentum technology for two-nucleon interactions 



In this section, we give an overview of the equations and techniques used to derive low-momentum NN 
interactions (with many-body interactions treated in Section H]). We use the equations as a guide to 
how renormahzation group methods decouple low- and high-energy degrees of freedom but refer to the 
literature for the more technical details. 



3.1 Sharp cutoff Mow/fc 

Imposing a sharp cutoff A in relative momentum is the most direct way to limit the resolution of an NN 
potential by excluding high-momentum modes. But to incorporate the relevant details from such modes 
into low-momentum interactions, they must be integrated out rather than simply truncated. The Viowfc 
approach [5|[58|l87] does this by demanding that the (half-on-shell) T matrix for an initial potential HvfN 
be unchanged in every NN partial wave as A is lowered, so that 

T{k ,k]k ) = V^^{k ,k) + -V p dp , (8) 

vr Jo k^ - p^ 

.A nj , 2^ f''V,l,{k',p)T{p,k;e^ 



= VLk{k\k) + -V ;r T ' ' rdp, (9) 

for all fc, k' < A. Note that this means that V^ow/c only has momentum components below the cutoff but 
the invariance of T{k', k; k"^) ensures that observables for these momenta (such as phase shifts and the 
deuteron binding energy) are preserved. 

Imposing dT{k', k; k'^)/dA = on Eq. gives the RG equation 

dy^ 2 V,t,{k',A)THA,k-A^) 



Integrating Eq. f lTOj) with the large cutoff initial conditioio = is equivalent to first re- 

summing high-momentum ladders in an energy-dependent effective interaction, which is the solution to 
the two-body Bloch- Horowitz equation in momentum space with the projector Q = ^ q^dq \q) and 
subsequently trading energy dependence for momentum dependence by using the equations of motion. 
This two-step procedure is also equivalent to the commonly-used Lee-Suzuki transformation method 
used to construct energy- independent effective interactions [51 I89149T] . Note, however, that the RG 
approach differs from Lee-Suzuki as the Q-space block of the effective Hamiltonian is set to zero. As a 
result of this equivalence (see Ref. [88]), the l^owfc interactions can be constructed in energy-independent 
form either using model-space methods (such as Lee-Suzuki or Okubo [02] transformations) or through 
the RG treatment. In practice, most calculations of Vfowfc have used the model-space methodfl] (followed 
by an Okubo Hermitization [5l[6]) rather than the differential equations of the RG because they are 
more robust numerically. This is true for both sharp and smooth cutoff versions of Mowfc- 

The evolution to lower cutoffs shifts contributions from the sum over intermediate states to the 
interactions, just as RG equations in quantum field theory shift strength from loop integrals to coupling 
constants. The change in the interactions include subleading contact interactions to all orders. Because 
only large momentum modes are decoupled, long-range pion exchanges are preserved. As documented 
in Section [273] different initial Vnn interactions collapse to the same nearly universal Viowfc at sufficiently 
small cutoffs (A ~ 2.0 fm"^). 



^The derivation of Eq. (jlOl) is slightly more complicated in the presence of bound states, although the final result 
remains the same [88] . 

^"The cutoff Aoo is chosen large enough so that the initial potential Vnn is effectively zero for k > Aoo- This can be 
very large for phenomenological potentials, for example, Aoo ^ 25 fm""'^ for the Argonne fig potential |18j . 
^^See Refs. [331121] for a three-dimensional formulation that avoids partial- wave decompositions. 
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3.2 Smooth regulators for V\o^ 



k 



The sharp cutoff RG equation for I^iowA; takes a particularly simple form that preserves two-body observ- 
ables for all momenta up to the cutoff. However, a sharp momentum cutoff leads to cusp-like behavior 



close to the cutoff in some channels and for the deuteron wave function (see inset in Fig. 25(a) ), which 
becomes increasingly evident as the cutoff is lowered below 2fm~^. In some applications, this leads to 
slow convergence, for example at the 10-100 keV level in few-body calculations using harmonic-oscillator 



bases, as shown in Fig. 26(b) One would expect that such calculations for the deuteron and the tri- 
ton, which are particularly low-energy bound states, should show improvement for cutoffs well below 
2fm~^, but instead a degradation was observed in Ref. [10]. This result was attributed to the use of 
sharp cutoffs, and subsequent studies [HllIT] have shown that these problems are alleviated by using a 
smooth-cutoff version of Viowk- 

The sharp cutoff l^owfc interactions were originally constructed using model-space methods (such as 
Lee-Suzuki [HEH] or Okubo [02] transformations) that rely on the definition of orthogonal projectors 
P and Q, such that P + Q = 1 and PQ = QP = 0. While smooth cutoffs seem incompatible with 
methods requiring PQ = 0, it is not a conceptual problem for the RG approach. Indeed, there is an 
appreciable literature on smooth-cutoff regulators for applications of the functional or exact RG [95] . 
The functional RG keeps invariant the full generating functional, which translates into preserving all 
matrix elements of the fully off-shell inter-nucleon T matrix. While this straightforwardly leads to RG 
equations, it also implies an energy- dependent interaction, which is undesirable for practical few- and 
many-body calculations. This conflict can be resolved by a generalization of the energy-independent 
RG equation derived in Section [3?T] to smooth cutoffs [6]. It is convenient to first define the partial-wave 
Viowk and the corresponding Tiowfc matrix in terms of a reduced potential v and a reduced t matrix as 

V,o^k{k',k) = f{k')v{k',k)f{k), (11) 
Tio„fc(fc', k; e) = f{k') t{k', k- P) f{k) . (12) 

where f{k) is a smooth cutoff function satisfying 

/(A;)^l and /(A;)^0. (13) 

Representative examples are exponentials f{k) = exp[— (/c^/A^)"°''p] and Fermi-Dirac functions f{k) = 
1/(1 -|- e*^'^^"^^^/'^^). The reduced half-on-shell t matrix obeys a Lippmann-Schwinger equation with loop 
integrals smoothly cut off by p (p) , 

t{k ,k]k^) = v{k ,k) + - I p^dp — . (14) 

TT Jo k^ - p^ 

Analogous to the RG derivation for a sharp cutoff, we impose that the reduced half-on-shell t matrix 
is independent of the cutoff, dt{k' , k; k'^)/dA = 0. This choice preserves the on-shell t matrix while also 
maintaining energy independence. The resulting RG equation is 



d 2 v{k\p)j^[p{pMp,k-y) 



which describes the evolution of the reduced low-momentum interaction with the cutoff. If one takes 
an energy-independent NN potential as the large-cutoff initial condition and numerically integrates 
Eq. f lT^ . then the resulting v preserves the half-on-shell T^n matrix for all external momenta, t(/c', fc; k"^) = 
T^^{k' , k] k"^) . Therefore, V\ovik{k\k) = f{k')v{k',k)f{k) preserves the low-momentum half-on-shell 
Tnn matrix up to factors of the smooth cutoff function. In the limit f{p) — )• 6{A — p), the RG equation 
for a sharp cutoff, Eq. (fTOl) . is recovered. 
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Figure 25: (a) S-state and (b) D-state components of the deuteron wave function in momentum space 
{u{k) and w{k) respectively) for the N'^LO potential of Ref. [20] and evolved Vfowik using smooth and 
sharp cutoffs at A = 2.0 fm~^ and 1.5 fm~^. 



The solution to Eq. f[T^ is non- Hermit ian. However, one can show [0] that a simple generalization 
of the Okubo Hermitization transformation to smooth cutoffs is obtained by symmetrizing the smooth- 
cutoff RG equation, Eq. f[T^ . to obtain 



d 1 /"°° 

— v{k',k) = - p^dp 
dA Tc Jq 



v{k',pWP{p)]t{p,k;p') ^ t{k',p;p')^[p{p)]v{p,k) 



(16) 



Here, the Hermitian low- momentum interaction V^iowA; preserves the low-momentum fully-on-shell Tnn 
matrix, up to factors of the regulator function Tio^k{k,k;k^) = p{k)Tj<ff^{k,k; k"^), and the deuteron 
binding energy. Note that demanding dt{k, k; k'^)/dA does not imply a "unique" Hermitian RG equation. 
In conventional model-space methods, this non-uniqueness corresponds to the choice in the Hermiti- 
zation procedure (or to residual unitary transformations in the P-space block). In RG language, this 
can be viewed as demanding dt{k', k; k'^)/dA = {k^ — k'^) $(fc', k), where $(A;', k) is any function satis- 
fying limfc^fc/(A;^ — k'"^) ^{k',k) = 0. The above RG equation derived from the Okubo transformation 
corresponds to a particular choice for ^{k', k). 

In practice, the numerical solution of Eq. f lTB]) is slowed by the t matrix evaluations involved at 
each step. In addition, the RG equation involves two-dimensional interpolations and principal-value 
integrals over narrowly peaked functions making it easy to introduce small errors at each step that 
can accumulate as the cutoff is lowered. Therefore, an alternative 3-step procedure to construct a 
low-momentum, energy-independent interaction with smooth cutoffs that is better suited for numerical 
calculations has been formulated (see Ref. [6] for details): 

1. Evolve a large-cutoff potential to a lower, smooth cutoff while preserving the full off-shell T 
matrix modulo factors of the smooth cutoff functions. This generates an energy- dependent low- 
momentum interaction, and amounts to solving the Bloch-Horowitz equation with a smooth cutoff. 

2. Convert the energy dependence to momentum dependence using equations of motion, which results 
in a non-Hermitian smooth-cutoff interaction that preserves the low-momentum half-on-shell T- 
matrix modulo factors of the smooth cutoff functions. 
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(a) (b) 

Figure 26: (a) The relative error in the deuteron binding energy Ed as a function of the size of the 
oscillator space for sharp cutoff and various smooth regulators, (b) The triton binding energy Et cal- 
culated from a direct diagonalization in a harmonic-oscillator basis of the low-momentum Hamiltonian 
evolved from the Argonne vi^ potential with cutoff A = 2 fm~ , as a function of the size of the oscillator 
space (A^max^ excitations). The open circles are calculated with a sharp cutoff for a fixed oscillator 
parameter h while the filled ones correspond to optimizing h at each A^max- The dashed line indicates the 
exact Faddeev result [S] , which shows the slow convergence at the 100 kev level of the diagonalization 
for the sharp cutoff. The squares are for a smooth Fermi-Dirac regulator and the triangles for a smooth 
exponential regulator, which each solve the convergence problem. For details see Ref. j6j. 



3. Perform a similarity transformation to Hermitize the low-momentum interaction. The result- 
ing low-momentum interaction is energy-independent and Hermitian, and preserves the low- 
momentum fully on-shell T-matrix up to factors of the smooth cutoff functions. 

It has been verified numerically that Eq. f fT6|) and the 3-step procedure give very similar results, and 
both tend towards the sharp cutoff (with Okubo Hermitization) calculated using the Lee-Suzuki 
method in the limit of a sharp cutoff. However, due to its computational advantages the 3-step method 
is used in most applications, although Eq. flTB]) has also been used in practical calculations [TTIISSIISZ]- 
The momentum space deuteron wave functions for smooth and sharp cutoffs are contrasted in Fig. [23 
We follow the notation of Ref. [71], with S- and D-state components denoted in coordinate space by 
u and w respectively, and with tildes in momentum space. The sharp-cutoff wave functions develop 



cusp-like structures in momentum space below 2fm~^ (see inset in Fig. 25(a)), which are removed by 
the Fermi-Dirac (or any other smooth) regulator. As with the sharp cutoff Mowfc; the smooth-cutoff 
evolution removes the short-range correlation "wound" in the coordinate-space deuteron wave function 
at lower cutoffs, and only modifies short- distance physics as indicated by the weak cutoff dependence of 



long-range operators like the "bare" quadrupole moment in the deuteron, see Fig. 16(b) These features 
lead to more perturbative behavior in nuclear matter |9] as well as in few-body systems [Tl][T2]. 

The primary motivation for smooth cutoffs was to remedy the slow convergence at the 10 keV level 
in the deuteron and at the 100 keV level in the triton, when calculated in a harmonic-oscillator basis. In 



Fig. 26(a) , the relative error in the binding energy of the deuteron (with respect to the converged result) 
is shown as a function of the size of the oscillator space (A^max ^ excitations) for a range of smooth 
regulators. The slow convergence is evident for the sharp cutoff, where the error is below the percent 
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level only for the largest space. The Fermi-Dirac regulator eliminates this problem, with improved 
convergence by increasing the smoothness parameter e from 0.2 fm~^ to 0.5 fm~^. Very similar results 
are obtained for the exponential regulator, or any other smooth regulator function. The convergence is 



also greatly improved for the triton. In Fig. 26(b) , the triton binding energy is plotted as a function of 
the size of the oscillator space. For efficiency, convergence for the smallest possible space is desirable, as 
the computational cost grows rapidly with Nraa,x and for larger systems. Convergence at the keV level 
is achieved by all exponential regulators with nexp ^ 4 (rzexp = 8 is shown) soon after A^max = 20. The 
consequence in moving from sharp to increasingly smooth regulators is also seen from the Fermi-Dirac 
regulators, where e = 0.5 fm~^ yields very satisfactory results. 

3.3 Similarity renormalization group (SRG) 

An alternative path to decoupling high-momentum from low-momentum physics is the similarity renor- 
malization group (SRG), which is based on a continuous sequence of unitary transformations that 
suppress off-diagonal matrix elements, driving the Hamiltonian towards a band-diagonal form [551 - 1102] . 
The SRG potentials are automatically energy independent and have the feature that high-energy phase 
shifts (and other high-energy NN observables) , while typically highly model dependent (see Fig. [1]), are 
preserved, unlike the case with Viowk as usually implemented. Most important, the same transforma- 
tions renormalize all operators, including many-body operators, and the class of transformations can 
be tailored for effectiveness in particular problems. 

In Ref . [7] , the first application of SRG transformations to NN interactions using the fiow equation 
formalism of Wegner [99 |I102] was made. The evolution or fiow of the Hamiltonian with a parameter s 
is a series of unitary transformations, 

H, = UMUl=T,,, + Vs, (17) 



where T^e\ is the relative kinetic energy and H = Trei -|- ^ is the initial Hamiltonian. Equation ([Tj 
defines the evolved potential Vg, with T^^i taken to be independent of s. Then Hs evolves according to 



^ = [Vs,H,], (18) 

with 

Vs = ^ Ul = -vl . (19) 

Choosing rjs specifies the transformation, which is taken as the commutator of an operator, Gg, with 
the Hamiltonian, 

r,s = [Gs,Hs], (20) 

so that 

[[Gs,H,],H,]. (21) 



ds 

Applications to nuclear forces have used Gg = Trd [7J, but one could also use momentum-diagonal 
operators such as Trei^, or the running diagonal Hamiltonian Hd, as advocated by Wegner [99] (see also 
Ref. [103] ). Taking the simplest choice Gs = Trei, the fiow equation in a partial wave takes the form 



dVs{k,k') 
ds 



9 1'°° 

-{k^ -k'^fVs{k,k') + - / q^dq{e + k'^-2q^)Vs{k,q)Vs{q,k'). (22) 

7^ Jo 

For matrix elements far from the diagonal, the first term on the right side of Eq. (1221) evidently dominates 
and exponentially suppresses these elements as s increases, 

V,{k,k')^Vs=o{k,k')e-<'''-'"^\ (23) 
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Figure 27: Contour plots of momentum-space matrix elements for the SRG evolution with A in the 
(a) ^So, (b) S-wave part of the ^Si-^Di and (c) ^Pi channels. The initial potential in (a) is the (A = 
600 Me V) N^LO potential [20] and in (b) and (c) the N^LO potential with A/A = 500/600 MeV 



It is convenient to switch to the flow variable X = s which has units of fm ^, because Eq. (1231) 
shows that A is a measure of the resulting diagonal width of Vg in momentum space. 

The evolution of the Hamiltonian according to Eq. fl22|) as s increases (or A decreases) is illustrated 
in Fig. [271 using two initial chiral EFT potentials. On top is ^Sq starting from the harder (A = 600 MeV) 
N^LO potential of Ref. [20], which has significant strength near the high-momentum diagonal, in the 
middle is the S-wave part of the ^Si-'^Di channel starting from one of the potentials of Ref. [H], which 
has more far off-diagonal strength initially and comparatively weaker higher-momentum strength on the 
diagonal, and on bottom is ^Pi with that same potential. Each of these examples show the characteristic 
features of the evolution in A, namely the systematic suppression of off-diagonal strength, as anticipated, 
with the width of the diagonal scaling as A^. 

The SRG-evolved interactions share key similarities (universality, increased perturbativeness, weaker 
correlations, etc.) with the smooth-cutoff l^owfc potentials, even though the decoupling of low and high 
momenta is achieved in a somewhat different manner. As A is lowered, different initial potentials flow 
to similar forms at low momentum (while remaining distinct at higher momentum), with the universal 
low-momentum parts numerically similar to the V^ow/c potentials. These observations are illustrated in 
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Figure 28: (a) Low-momentum universality for momentum-space matrix elements of the evolved SRG 
potentials at A = 2fm~^ for ^Sq (top, diagonal elements) and ^Si (bottom, off-diagonal elements). Also 
shown is the Vfowfc potential for a smooth regulator with A = 2fm~^ and ngxp = 4. (b) Largest repulsive 
Weinberg eigenvalues in the ^Sq and ^Si-^Di channels as a function of A, with initial potentials as in 
Fig. EZl For details see Ref. [7]. 
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Figure 29: 
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a) Absolute error versus A of the deuteron binding energy from a variational calculation in 
basis of harmonic oscillators (A^max^ excitations), (b) Variational triton binding energy 
A with NN interactions only, as a function of A^max- For details see Ref. [7]. 



Fig. 28(a) for the diagonal and off-diagonal matrix elements. 

Similarly, the nonperturbative features associated with strong short-range repulsion and strong 
short-range tensor forces are substantially softened as we evolve to lower A. The largest repulsive 
Weinberg eigenvalues (see Section l27il) for E = are shown as a function of A in Fig. 28(b) for the ^Sq 
and ^Si-^Di channels. In both channels, the large eigenvalues decrease rapidly as A evolves to 2fm~^ 
and below, as observed with the Viowfc evolution in Section [2^ The more perturbative potentials at 
lower A induce weaker short-range correlations in few- and many-body wave functions, which leads to 
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(a) (b) 

Figure 30: (a) and (b): Comparison of momentum-space Viowfc (left) and SRG block-diagonal (right) 
potentials with A = 2fm~^ evolved from the N^LO ^Si potential of Ref. [20]. The color and z axes are 
in fm. For details see Ref. |104] . 



greatly improved convergence in variational calculations. This is illustrated in Figs. 29(a) and 29(b) via 
calculations of the deuteron and triton binding energies by diagonalization in a harmonic-oscillator basis. 
The improvement in convergence is similar to that found with smoothly-cutoff Vio„fc interactions [6]. 

In Fig. 29(b) , the calculations for different A converge to different values for the triton binding 
energy. This reflects the contributions of the omitted (and evolving) three-body interactions, and 
follows a similar pattern to that seen with NN-only Viow/t calculations [6l|8]. The consistent evolution 
of many-body forces is an important issue for low-momentum interactions. For the SRG, the evolution 
of 3N forces is readily practical, as discussed in Section 14.21 



3.4 Mow/c from SRG flow equations 

A powerful feature of the SRG is that the generator Gs can be tailored to decouple high- and low- 
momentum physics in different ways |104j . Block-diagonal decoupling of the sharp Mowfc form can be 
generated using SRG flow equations by choosing a block-diagonal flow operator |105[I106] . 











(24) 



with projection operators P and Q = \ — P . In a partial-wave momentum representation, P and Q are 
step functions defined by a sharp cutoff A on relative momenta. This choice for G^, which means that 
r\s is non-zero only where is zero, suppresses off-diagonal matrix elements such that the Hamiltonian 
approaches a block-diagonal form as s increases. If one considers a measure of the off-diagonal coupling 
of the Hamiltonian, 

Ti\{QE,P)\QE,P)\ = Ti[PHsQHsP] ^ , (25) 
then its derivative is easily evaluated by applying the SRG equation, Eq. fl2T]) : 



^ Tr[PHsQHsP] = Ti[Pr]sQ{QH,QHsP - QH.PHsP)] + Ti[{PH,PH,Q - PHsQHsQ)Qr]sP] 



= -2Ti[iQrisPyiQr]sP)]^0. 
Thus, the off-diagonal QHgP block will decrease in general as s increases [10511106] . 



(26) 



The right panels of Figs. 30(a) and 30(b) result from evolving the N^LO potential of Ref. pO] using 
the block- diagonal Gg of Eq. (!24|) with A = 2fm~^ to A = 0.5 fm~^. The agreement between Vio^k and 
SRG potentials for momenta below A is striking, where a similar degree of universality is found in other 
channels. Deriving an explicit connection between these approaches is the topic of ongoing research. 



32 



^0 2 4 6 8 10 




2 4 6 8 10 ^0 2 4 6 8 10 



10 



10 



D 2 4 


6 8 10 









Figure 31: Evolution of the ^Si partial wave using the SRG block-diagonal flow equation with A = 2 fm ^ 
at A = 4,3,2, and 1 fm"^ [lOl]. The initial N^LO potential is from Ref. [20]. The axes are in units of 



k from — llfm and the color scale ranges from —0.5 to +0.5 fm as in Fig. 30(a 



The evolution with A of the ^Si channel is shown in Fig. |3T] (for results in the ^Pi channel see 
Ref. |104] ). The evolution of the "off- diagonal" matrix elements (meaning those outside the PHgP and 
QHsQ blocks) can be roughly understood from the dominance of the kinetic energy on the diagonal. 
Let the indices p and q run over indices of the momentum states in the P- and Q-spaces, respectively. 
To good approximation we can replace PHgP and QHgQ by their eigenvalues Ep and Eq in the SRG 
equations, yielding |105[I106] 



ripq ~ ^pq ^q ^p Vpq 



-{Ep - Eq) r]pq and 7]pq ^ Ep hpq - hpq Eq = {Ep 

Combining these two results, we have the evolution of any off-diagonal matrix element: 

-{Ep-E,?h 



P^q) hpq 



ds 



(27) 



(2J 



For two-nucleon interactions we can replace the eigenvalues by those for the relative kinetic energy, 
which gives the explicit solution 

hpq{s)^hpq{0)e-'^p"-'^"K (29) 

Thus the off-diagonal matrix elements go to zero with the energy differences just like with the SRG 
with Trei; one can see the width of order I/a/s = A^ in the k'^ plot of the evolving potential in Fig. EH 
(this is even clearer in other partial waves, see Ref. |104j ). 

While in principle the evolution to a sharp block-diagonal form means going to s = oo (A = 0), 
in practice one must only take s as large as needed to quantitatively achieve the decoupling implied 
by Eq. ( l29|) . Furthermore, this should hold for more general definitions of P and Q. To smooth out 
the cutoff, one can introduce a smooth regulator /a, which can be taken as an exponential /a (A;) = 



e ('^Va^ _ Pqj, Viowk interactions, typical values used are 



n, 



exp 



4, 6, and 8 (the latter is considerably 



sharper but still numerically robust). By replacing the generator with 

Gs = fAHjA + {l-h)Hs{l-fA), 

one obtains a smooth block- diagonal potential (for examples, see Ref. |104j ). 



(30) 



3.5 Related methods 

There are various other methods used to transform nuclear Hamiltonians for use in few- and many-body 
calculations. Although they are not cast in the RG language, many of the same principles and features 
discussed in Section [2] apply. While space does not permit detailed discussions of each, here we give a 
brief synopsis and comparison, with pointers to the literature for more information. 
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(a) Fav18 (b) VsRG (c) 1'uCOm(SRG gen.) (d) X'ucOM(var. opt.) 




Figure 32: Comparison of momentum-space matrix elements of ^Sq potentials (in units of MeVfm'^) 
for (a) Argonne fig, (b) SRG-evolved with A = 2.4 fm~^, (c) UCOM-transformed using SRG-generated 
correlators (A = 2.25 fm~^), and (d) UCOM-transformed using variationally optimized correlators |107j . 
Reprinted with permission from R. Roth et al. |107j . copyright (2008) by the American Physical Society. 

Bloch-Horowitz. The energy-dependent Bloch-Horowitz (BH) equation is used as an intermedi- 
ate step in deriving energy- independent sharp and smooth-cutoff Viowfc interactions. The BH method 
has also been used directly in other approaches, where the energy is determined self-consistently. The 
harmonic-oscillator-based effective theory |108H111] is one such example. While most detailed appli- 
cations have been limited to A = 2, 3 because energy-dependent interactions are difficult to handle in 
other approaches, perturbative expansions and clever resummations may overcome this hurdle. 

UCOM. Feldmeier, Neff, Roth and collaborators have developed a method using unitary trans- 
formations designed to remove short-range "hard core" and tensor correlations |101llll2] - lll4j . The 
approach is called the Unitary Correlation Operator Method (UCOM). The result is a soft Hamilto- 
nian that shares favorable features of the RG-based interactions. They have also studied parallels (and 
some distinctions) of UCOM to the SRC approach [TOTlinS], as shown in Fig. [31 Using a UCOM NN 
potential, which is phase equivalent to a given initial Hamiltonian, they have investigated ground-state 
energies and radii of many nuclei (including heavy nuclei such as ^^'^Pb) in Hartree-Fock plus many-body 
perturbation theory. Preliminary results from supplementing the NN potential by a short-range three- 
body force with fitted strength are encouraging |116] . As originally formulated, the UCOM procedure 
is difficult to carry out beyond the two-body level. However, the connections to the SRG approach 
established in |107[lll5j may help to transform 3N forces. 

NCSM Lee-Suzuki. The Lee-Suzuki transformations used in the No-Core-Shell-Model (NCSM) 
are constructed to decouple low from high energy within a harmonic-oscillator model space [TIT] . A 
cluster expansion is used to organize and truncate the many-body components. This approach has 
been successfully used to describe nuclear spectra, especially in p-shell nuclei starting from chiral EFT 
interactions [57 |I118] . A detailed overview of this approach including recent developments is given in 
Ref. [nS]. As usually constructed, the Lee-Suzuki approach is not variational and does not always show 
smooth convergence with respect to the model-space size. Unlike the momentum-space RG evolution, 
it is not clear if evolving in the harmonic-oscillator basis is problematic for long-range operators, as 
indicated by the sizable effective three- and higher-body contributions to B{E2) values [65] . 

UMOA. In the Unitary Model Operator Approach (UMOA), a Hermitian, energy- independent 
effective interaction is derived by Okubo-type unitary transformations of an initial Hamiltonian (see 
Ref. |119j and references therein). The transformation is applied twice to make the effective interaction 
more suitable to larger nuclei. In the first transformation, the short-range repulsion is softened by 
partitioning the P- and Q-space with an energy cutoff on two-body harmonic-oscillator states. In the 
second transformation, the softened interaction is transformed so that there are no matrix elements 
between the unperturbed ground state and two-particle-two-hole states for a given closed-shell nucleus. 
Calculations have been performed for stable nuclei around ^^O [120] . A-hypernuclei |121] and for ground- 
state and single-particle energies of ^°Ca and ^^Ni |122] . 



34 



k(fm"'') k(fm"^) k(fm"^) k(fm"^) 




Figure 33: G matrix at saturation density for the Argonne fis potential [18] (left panels) and the N^LO 
potential of Ref. [20] (right panels) in the ^Si channel. Each set of four panels are a) initial potential, 
b) potential evolved by the SRG to A = 2fm~^, c) G matrix based on a), and d) G matrix based on b). 

G matrix. The G matrix is the in-medium sum of ladder diagrams to form a two-body effective 
interaction from a starting NN potential. This formalism has been used in calculations of nuclear mat- 
ter |123j . for effective interactions in nuclei [91], and in self-consistent Green's function approaches |124] . 
It has been suggested that V^owfc NN interactions are essentially the same as a G matrix but there are 
substantial differences. First, Mowfc at a given cutoff A is a free-space interaction, equivalent in practice 
to any low-energy Hamiltonian (with no "starting energy" dependence). Second, significant coupling 
between low and high momenta remains in the G matrix. This is demonstrated in Fig. |33l which shows 
the Argonne fig and an N^LO potential in the ^Si channel compared to low-momentum interactions 
(SRG evolved to A = 2fm~^, similar results are found for Viowfc interactions) and to the G matrices 
obtained from each. (The qualitative conclusions do not depend on the particular choice of starting en- 
ergy.) We observe that the summation into the G matrix has a relatively small effect on low-momentum 
interactions, in stark contrast to the Argonne fig potential case. The resulting low-momentum matrix 
elements are indeed similar, but it is also apparent that there is substantial off-diagonal strength in 
the G matrix in this case. This coupling to high momenta prevents the expansion of nuclear matter 
properties from being perturbative in the G matrix, necessitating a non-perturbative scheme such as 
the hole-line expansion. One might expect the softer N^LO ^Si potential to avoid this problem, but we 
see that the off-diagonal strength, while limited to lower momentum, is actually increased. 

4 Many-body interactions and operators 

When evolving nuclear interactions to lower resolution, it is inevitable that many-body interactions 
and operators are induced even if initially absent. This might be cause for alarm if nuclei could be 
accurately calculated with interactions truncated at the two-body level, as was assumed for part of the 
history of nuclear structure calculations. However, chiral EFT reveals the natural scale and hierarchy of 
many-body forces, which dictates their inclusion in modern calculations of nuclei and nucleonic matter. 
Thus the real concern is whether this hierarchy is maintained as nuclear interactions are evolved. In 
this section, we review the current status of RG technology to include many-body interactions and 
operators and the currently known impact on the hierarchy. 
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4.1 Three- nucleon interactions 



Three-nucleon interactions are a frontier. They are crucial for binding energies and radii, they play 
a central role for spin-orbit effects, spin dependencies, for few-body scattering and the evolution of 
nuclear structure with isospin, and they drive the density dependence of nucleonic matter (see Sections [5] 
and[6j) [73]. Three-nucleon interactions are also required for renormalization |125pi26j . The construction 
of 3N forces based on chiral EFT provides a systematic organization of the physics and an operator 
basis that can be used to approximate the evolution of low-momentum 3N interactions. 



In chiral EFT without explicit A isobars, 3N forces first enter at N^LO (see Fig. 4(a) ) and contain a 
long-range 27r-exchange part Vc, an intermediate-range Ivr-exchange part Vd and a short-range contact 
interaction Ve [T271 fT28]: 



(31) 
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The 27r-exchange interaction is given by 





2 V2/. 



(32) 



where = k ■ — kj denotes the difference of initial and final nucleon momenta and k = 1,2, 3) and 



4:Ciml 2c3 

J IT J TV 



(33) 



while the Ivr-exchange and contact interactions are given respectively by 



Vd 
Ve 



8/' /'Ax + < 



[Ti ■ Tj) [CTi ■ q 



■J J ' 



ce 



ii^k 



(34) 
(35) 



Typical values for applying Eqs. ([32])-([35]) are qa = 1-29, U = 92.4 MeV, = 138.04 MeV and 
= 700 MeV. In the RG calculations based on chiral EFT interactions discussed here, the 3N force 
contributions are regulated as in Ref. [i6j using 



/r(p, q) = exp 



(p^ + 3gV4)- 



(36) 



with a 3N cutoff A3NF that is allowed to vary independently of the NN cutoff. Here, p and q are Jacobi 
momenta. The exchange terms of the 3N force are included by means of the antisymmetrizer 



A 



123 



{1 + P12P23 + Pl3P23){l - P'. 



23) 



13 



-P23 + P12P23 



P3P 



13-' 23 



(37) 



where Pij is the exchange operator for spin, isospin and momenta of nucleons i and j. The regulator 
/r(P; q) is totally symmetric when expressed in the nucleon momenta kj, and thus the direct and 
exchange terms contain the same regulator. 

The low-energy q couplings for the long-range parts relate vrN, NN and 3N interactions, which means 
they can be fitted in different processes. The determination of the Cj couplings from ttN scattering agrees 
reasonably well with the extraction from NN partial waves ^1] . However, there are large uncertainties in 
the values (and some controversy, see Ref. |129] ). More accurate determinations are needed for nuclear 
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1.0 


21.06 


-28.62 


0.02 


0.11 


-1.06 


38.11 


-62.18 


0.10 


0.54 


-4.87 


0.08 


0.55 


1.3 


25.71 


-34.14 


0.01 


1.39 


-1.46 


50.14 


-78.86 


0.19 


8.08 


-7.83 


0.10 


0.63 


1.6 


28.45 


-37.04 


-0.11 


0.55 


-0.32 


57.01 


-86.82 


-0.14 


3.61 


-1.94 


0.04 


0.67 


1.9 


30.25 


-38.66 


-0.48 


-0.50 


0.90 


60.84 


-89.50 


-1.83 


-3.48 


5.68 


0.06 


0.70 


2.5(a) 


33.30 


-40.94 


-2.22 


-0.11 


1.49 


67.56 


-90.97 


-11.06 


-0.41 


6.62 


0.12 


0.74 


2.5(6) 


33.51 


-41.29 


-2.26 


-1.42 


2.97 


68.03 


-92.86 


-11.22 


-8.67 


16.45 


0.18 


0.74 


3.0(*) 


36.98 


-43.91 


-4.49 


-0.73 


3.67 


78.77 


-99.03 


-22.82 


-2.63 


16.95 


0.23 


0.80 



Table 1: Expectation values of the kinetic energy T, V^owfc and the different 3N contributions (long-range 
27r-exchange part V^, Ivr-exchange part Yd and contact interaction Ve) for and ^He. All energies are 
in MeV and momenta are in fm~^. (a) and (6) denote two possible solutions for A = 2.5 fm~ and (*) 
indicates that the ^He fit is approximate, for details see Ref. [8]. Also listed is the ratio of maximum 
3N to V^owfc contribution and an average relative momentum fcrms [HS]- 
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Figure 34: The and ^He radii are approximately cutoff independent with low-momentum NN and 
3N interactions [73lll33j . 
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structure, where C3 and C4 are particularly important (see Fig. 44(b) for an example of the impact 
of the Ci uncertainties). The remaining cd and ce couplings are usually fit to the binding energy 
and another observable in A ^ 3. These leading 3N interactions generally improve the agreement of 
theory with experiment in nucleon-deuteron scattering [311 136] and in the spectra of light nuclei [57] 
(see Figs. 5(a) and[H]). The subleading 3N interactions at N^LO are parameter free and several parts 
have been calculated to date |13Uyi31j . The subleading 3N interactions at N^LO are parameter free and 
several parts have been calculated to date |130[ll31j . including two-pion-one-pion-exchange and pion- 
ring diagrams, spin-orbit forces, and contributions that involve NN contacts [31]. One consequence is 
that operator structures of the form in Eq. ( 132|) are generated at N'^LO, effectively weakening the values 
of the Ci couplings [31] . 

Instead of explicitly evolving 3N interactions to low momentum, one can take the N^LO 3N forces 
as a truncated basis and assume that the Cj couplings are not modified by the RG evolution. This is 
the strategy adopted in Refs. [8l|9l[T6]. In this scheme, V3n(A) is constructed by fitting the cd and ce 
couplings to the '^H binding energy and the "^He binding energy [8] or matter radius for a range of 
cutoff values. When the ^He radius is fit, the predicted ^He binding energies are very reasonable. Fit 
values for ce, and ce are found to be natural and are tabulated in Refs. [S] (with Argonne fis [IS] as 
initial NN potential) and [IS] (with the N^LO potential of Ref. [2D])- For cutoffs A < 2.0fm"\ there 
are linear dependencies in the fitting, which are consistent with a perturbative Vd and Ve contribution: 
E^s = -^(Viowfe + Vc) + cei{D— term) + C£;(£'— term) (where (. . .) denotes the matrix elements of the 
operators). This has been verified explicitly and also for Vc. 
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Thus, when evolved to lower cutoffs, low-momentum 3N interactions become perturbative in light 
nuclei [Sl|T6j which means (\E''^^^|V3n|^^^'') ~ (^'•^^IV'snI^I''-^-'), where |vl/(")) are exact solutions including 
up to n-body forces. It is instructive to look at the expectation values of different contributions to the 
and ^He ground-state energies, as listed in Table [1] The expectation values are consistent with chiral 
EFT power counting estimates, which means that the hierarchy of many-body forces is maintained. 
Note also that the sign of the q contribution is cutoff dependent. This contribution is repulsive in 
nuclear and neutron matter for low-momentum interactions (see Sections 15. II and l5.2p . which is counter 
to the intuition from conventional NN and 3N potentials [5311134] . 

As discussed in Section flM (see Fig. neglecting 3N interactions leads to a universal correlation 
(empirically known as Tjon-line) between the and ^He binding energies [StlSl]. This is an example of 
how cutoff variation can provide lower bounds for theoretical uncertainties due to neglected many-body 
interactions or an incomplete many-body treatment. Another example is shown in Fig. [311 where the 
radii for '^H and ^He are shown as a function of Viowfc cutoff A. With 3N forces included, there is only a 
small cutoff dependence that shows the scale of the residual dependence on higher-order shorter-range 
many-body interactions. 



4.2 Three-nucleon force evolution 

While the use of a chiral EFT basis to approximate low-momentum 3N interactions is convenient and 
computationally efficient, it has not been verified against explicitly evolved three-body forces. This 
will be possible using SRG flow equations, which offer a tractable approach to evolving 3N forces (that 
avoids solving for any T matrices). The SRG flow equation, 

^ = [[G„H,],Hs], (38) 

can be applied directly in the three-particle space. The right side involves only the Hamiltonian and the 
generator Gg, which can be evaluated in a basis without solving bound state or scattering equations. 
A potential issue is the role of spectator nucleons; we consider two solutions: first a decoupling of the 
3N part in momentum-space representation and then a direct solution in a harmonic-oscillator basis. 
In both cases, we take Gg = Trei and return at the end to other choices. 

To show the basic idea in momentum space, we adopt a notation in which Vu means the two-body 
interaction between particles 1 and 2 while V123 is the irreducible three-body potential. We start with 
the Hamiltonian including up to three-body interactions (keeping in mind that higher-body interactions 
will be induced as we evolve in s but will not contribute to three-body systems): 

Hg = T,el + 1^12 + ^13 + V23 + Vu3 = T,el + K ■ (39) 

(Note: all of the potentials depend implicitly on s.) The relative kinetic energy operator Tj-ei can be 
decomposed in three ways: 

T,ei = Ti2 + Ts = Tis + T2 = T2-i + Ti , (40) 



and Tj commutes with 



[Ts, V^i2] = [^2, = [Ti, V23] = , (41) 



so the commutators of Trei with Vjk become [Tj-ei, V12] = [T12, V12] and similarly for V13 and V23. 

Because we define T^^i to be independent of s, the SRG flow equation, fl55]) . for the three-body 
Hamiltonian Hs simplifies to 

dVs _ dVu dV,, dV,s dV,,, _ 
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Figure 35: Ground-state energy of (a) '^H and (b) ^He as a function of A. For details see Ref. [85] 



with Vg defined by Eq. (1391) . The corresponding equations for each of the two-body potentials (which 
are completely determined by their evolved matrix elements in the two-body systems) are 



12 



ds 



[[Ti2, 1^12], + 1^12)], 



(43) 



and similarly for V13 and V23. After expanding Eq. (142|) using Eq. ( 139|) and the different decompositions 
of Trei, it is straightforward to show that the derivatives of two-body potentials on the left side cancel 
precisely with terms on the right side, leaving 



dV, 



123 



ds 



[[T12, ^12], (T3 + Fi3 + V^23 + ^123)] + [[Tl3, 1^13], (T2 + Vx2 + ^23 + ^123)] 
+ [[T23, ^23], (Ti + ^12 + Vx^ + 1^123)] + [[Trel, 1^123], . 



(44) 



The importance of these cancellations is that they eliminate the "dangerous" delta functions, which 
make setting up the integral equations for the three-body system problematic |135j . We emphasize that 
the s-dependence of the two-body potentials on the right side of Eq. ( 1441) is completely determined by 
solving the two-body problem in Eq. ( 122|) . This is in contrast to RG methods that run a cutoff on the 
total energy of the basis states (as in the Bloch-Horowitz or Lee-Suzuki approaches). Such methods 
generate "multi- valued" two-body interactions, in the sense that the RG evolution of two-body operators 
in A > 2 systems depends on the excitation energies of the unlinked spectator particles |109[ll36j . 

Further simplifications of Eq. fl441) follow from antisymmetrization and applying the Jacobi identity, 
but this form is sufficient to make clear that there are no disconnected pieces. The problem is thus 
reduced to the technical implementation of a momentum-space decomposition analogous to Eq. fl22|) . 
A diagrammatic approach is introduced in Refs. [137(1138] to handle this decomposition. Work is in 
progress on evolving 3N forces in momentum space. It has been verified that this formalism leaves 
eigenvalues invariant for three-particle systems described by simple model Hamiltonians, such as a 
two- level system of bosons |137] . 

To summarize, because only the Hamiltonian enters the SRG flow equations, there are no difficulties 
from having to solve T matrices (bound state plus scattering wave functions) in all three-body (including 
breakup) channels, as required by the analogous three-body Viowk evolution equations. In a momentum 
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(a) (b) 

Figure 36: (a) Ground-state energy of '^He as a function of the basis size A^max based on the N^LO 
potential of Ref. [20] and N^LO 3N forces from Ref. [139] . Unevolved (bare) results are compared to 
those obtained with a Lee-Suzuki (LS) NCSM transformation and based on SRG-evolved interactions 
to A = 2.0 fm~ at Ml = 28 and 36MeV. (b) Binding energy of ^He versus ^H. The Tjon line from 
phenomenological NN potentials (dotted) is compared with the trajectory of SRG energies when only NN 
interactions are kept (circles). When 3N interactions (initial and induced) are included, the trajectory 
lies close to experiment (cross) for A ^ 1.8 fm~^ (see inset). For details see Ref. [53] . 

basis the presence of spectator nucleons requires solving separate equations for each set of (v!^^^) matrix 
elements. But while it is natural to solve Eq. fl38|) in momentum space, it is an operator equation so 
one can use any convenient representation. A convenient choice is to evolve in a discrete basis, where 
spectators are handled without a decomposition. This results in coupled first-order differential equations 
for the matrix elements of the flowing Hamiltonian Hs where the right side of Eq. ( l38i) is evaluated 
using simple matrix multiplications. 

Such calculations have been recently performed in the Jacobi coordinate harmonic-oscillator basis of 
the No-Core Shell Model (NCSM) [65] . This is a translationally invariant, antisymmetric basis for each 
A, with a complete set of states up to a maximum excitation of N^^^hfl above the minimum energy 
configuration, where Q is the harmonic oscillator parameter. The 3N evolution used builds directly 
on Ref. [138] . which presents a one- dimensional implementation of the approach along with a general 
analysis of the evolving many-body hierarchy. 

One starts by evolving Hx in the A = 2 subsystem, which completely fixes the two-body matrix 
elements (V^^^V Next, by evolving Hx in the A = 3 subsystem, the combined two-plus-three-body 
matrix elements are determined. The three-body matrix elements are isolated by subtracting the 
evolved (V^ ) elements in the A = 3 basis [138] . Having obtained the separate NN and 3N matrix 
elements, they can be applied unchanged to any nucleus. The inclusion of 3N forces in the initial 
Hamiltonian is straightforward. If applied to A ^ 4 systems, four-body (and higher) forces will not be 
included and so the SRG transformation will be only approximately unitary. This allows one to study 
whether the decreasing hierarchy of many-body forces is maintained and whether the induced four-body 
contribution is unnaturally large. 

Results from Ref. [65] for the SRG evolution are shown in Fig. |35l The initial NN interaction is 
the (A = 500 MeV) N^LO potential of Ref. [20] while the initial 3N interaction is the N^LO 3N force 
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(discussed in Section I4.ip in the local form of Ref. |140] with cd and ce couplings fit to the average 
of the '^H and '^He binding energies and to triton beta decay according to Ref. |139] . Similar results 
are expected starting from other initial potentials because of the low-momentum universality. NCSM 
calculations with these initial interactions and the parameter set in Table I of Ref. |139j yield energies of 
-8.473(4) MeV for and -28.50(2) MeV for ^He compared with -8.482 MeV and -28.296 MeV from 
experiment, respectively. So there is a 20keV uncertainty in the calculation of ^He from incomplete 
convergence and a 200 keV discrepancy with experiment. The latter is consistent with the omission of 
three- and four-body interactions at N^LO. These provide a scale for assessing whether induced four- 
body contributions are important compared to other uncertainties. Figure [35] shows the SRG evolution 
as a function of A of the and ^He ground-state energies. When the induced three-body part is 
included, the ground-state energy is completely independent of A. For ^He, the deviation from A 
independence due to induced four-body forces is less than about 50keV for A ^ 1.5 fm~^. 



As seen in Fig. 36(a), the NCSM calculations with SRG-evolved interactions are variational and 
converge smoothly and rapidly from above with or without an initial three-body force. This shows the 
dramatic improvement in convergence rate compared to the initial chiral EFT interactions (circles), 
although the latter are relatively soft. Thus, once evolved, a much smaller A^max basis is adequate for 
a desired accuracy and extrapolating in Amax is also feasible (in contrast to the Lee-Suzuki NCSM 
results). The impact of evolving 3N forces is illustrated in Fig. 36(b) , where the binding energy of "^He 
is plotted against the binding energy of ^H. The experimental values define a point in this plane. The 
SRG NN-only results follow a trajectory that is analogous to the Tjon line (see Fig. |2^ . In contrast, the 
short SRG trajectory including 3N interactions highlights the small variations from omitted short-range 
three- and four-body forceso 

These results demonstrate that the SRG is a practical method to evolve 3N (and higher-body) forces 
in a harmonic-oscillator basis. Calculations of A ^ 4 nuclei including 3N forces show the same favorable 
convergence properties observed with low-momentum NN interactions, with a net induced four-body 
contribution m A = A that is smaller than the truncation errors of the initial chiral interactions. 
However, because of the stronger density dependence of four-nucleon forces, it will be important to 
monitor the size of the induced four-body contributions for heavier nuclei and nuclear matter. The SRG 
evolution is an alternative to the Lee-Suzuki NCSM transformation and the SRG-evolved harmonic- 
oscillator matrix elements can also be used with other many-body methods. The success of the block- 
diagonal SRG (see Section [3^ in evolving a V\o^k two-body interaction using flow equations motivates 
an extension to three- and higher-body forces (and operators). Work on this is in progress. 



4.3 In-medium SRG 

The SRG methods described so far carry out the evolution to low momentum at zero density in free 
space. A convenient feature of the free-space evolution is that it does not have to be repeated for each 
different nucleus or nuclear matter density; the evolved NN interaction is completely determined by 
solving the two-body SRG equation in free space, the evolved 3N interaction is completely determined 
by solving the three-body SRG equation in free space, and so on. 

An interesting alternative is to perform the SRG evolution directly in the A-body system of interest. 
In contrast to the free-space SRG, the in-medium evolution must be repeated for each nucleus or density. 
As with all RG methods, evolving to low momentum generates many-body interactions that must be 
truncated for practical calculations. However, unlike the free-space evolution, the in-medium SRG has 
the appealing feature that one can approximately evolve 3, ...,y4-body operators using only two-body 
machinery. The key to this simplification is the use of normal-ordering with respect to a finite-density 
reference state. That is, starting from the second-quantized Hamiltonian with two- and three-body 

^^Note that the trajectory for NN plus 3N-induced results would be a similarly small line at the N'^LO NN-only point. 
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interactions, 



H = J2Ti2a\a2 + (12| V|34)44a4a3 + 7^ 5Z (123|\/(3) |45g^^t^t^t^^^^^^ ^ ^45^ 

12 ^ 1234 ^ 123456 

all operators are normal-ordered with respect to a finite-density Fermi vacuum |$) (for example, the 
Hartree-Fock ground state or the non-interacting Fermi sea in nuclear matter), as opposed to the zero- 
particle vacuum. Wick's theorem can then be used to exactly write H as 

i7 = Eo+5^/i2{aIa2} + -^5^(12|r|34){al4a4a3} + ^ (123|r(3)|456){al44a6a504} , (46) 

12 ^ 1234 ^ 123456 

where the zero-, one-, and two-body normal-ordered terms are given by 

'"■i'^2 + ^ (123|V^^^|123)nin2n3 , (47) 

1 12 ■ 123 



/12 = T12 + 5^(lz|\/|2^)n, + i J2^Uj\W\2zj)n,nj , (48) 

i ij 

(12|r|34) = {12\V\3A) + ^(12z|1/(3)|34i)ni , (49) 



where Ui = 6{eF — Ei) denotes the sharp occupation numbers in the reference state, with Fermi level 
or Fermi energy ep. By construction, the normal-ordered strings of creation and annihilation operators 
obey ($|{a| ■ ■ ■ a„}|$) = 0. It is evident from Eqs. fl47l) - fl49l) that the coefficients of the normal-ordered 
zero-, one-, and two-body terms include contributions from the three-body interaction V"*^'^-' through 
sums over the occupied single-particle states in the reference state |$). Therefore, truncating the 
in-medium SRG equations to two-body normal- ordered operators will (approximately) evolve induced 
three- and higher-body interactions through the density-dependent coefficients of the zero-, one-, and 
two-body operators in Eq. ( H6l) . 

The in-medium SRG flow equations at the normal-ordered two-body level are obtained by evaluating 
dH/ds = [rj, H] with the normal-ordered Hamiltonian H = Eq+J+T and the SRG generator rj = ■q^^J^rj^^ 
(with one- and two-body terms) and neglecting three- and higher-body normal- ordered terms. For 
infinite matter, a natural generator choice is 77 = [/, F] in analogy with the free-space SRG. In this case, 
the explicit form of the SRG equations simplifies because 77^^ = and fij = fiSij. This leads to 

^ = ^X^(/l2 - /34)|ri234p^l"-2^3^4, (50) 



1234 



^ = (/4I - /23)|r4123p(^2'^3'^4 + niUsni) , (51) 



234 

'^ri234 



= -(/l2 - /34)^ ri234 + ^ X] (/l2 + /34 " 2/ab)ri2abra,,34(l - -«;,)+ ^ (ra„ - Tife) 

ab ab 

X |ral63r62a4 [{fal — fbs) " ifb2 — fai)] — ra263r61a4 [{fa2 — fbs) " {fbl — fai)] | , (52) 



where the single-particle indices refer to momentum states and include spin and isospin labels. While 
the in-medium SRG equations are of second order in the interactions, the fiow equations build up non- 
perturbative physics via the successive interference between the particle-particle and the two particle- 
hole channels in the SRG equation for F, Eq. ( 152|) . and between the two-particle-one-hole and two- 
hole-one-particle channels for /, Eq. ( [5T]) . In terms of diagrams, one can imagine iterating the SRG 
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A, [fm '] X [fm"'] 

(a) (b) 

Figure 37: In-medium SRG evolution of the energy per particle of (a) symmetric nuclear matter at 
a Fermi momentum kp = 1.4 fm~ and (b) neutron matter at kp = 1.35 fm~ |141] . The solid line 
denotes the Eq flow of the Hartree-Fock energy. The dashed line is the energy calculated in the ladder 
or Brueckner-Hartree-Fock approximation (BHF) using interactions evolved with the in-medium SRG 
at the two-body level, and the dotted line is the same many-body calculation based on NN-only SRG 
interactions evolved in free space. The initial NN interaction is the N^LO potential of Ref. [20]. 



equations in increments of 6s. At each additional increment 6s, the interactions from the previous step 
are inserted back into the right side of the SRG equations. Iterating this procedure, one sees that the 
SRG accumulates complicated particle-particle and particle-hole correlations to all orders. 

With this choice of generator rj = [/, F], the Hamiltonian is driven towards the diagonal. This 
means that Hartree-Fock becomes increasingly dominant with the off-diagonal F matrix elements being 
driven to zero. As with the free-space SRG, it is convenient for momentum-space evolution to switch 
to the flow variable A = s~^^^, which is a measure of the resulting band-diagonal width of F. In 
the limit A — )■ 0, Hartree-Fock becomes exact for the evolved Hamiltonian; the zero-body term, Eq, 
approaches the interacting ground-state energy, / approaches fully dressed single-particle energies, and 
the remaining diagonal matrix elements of F approach a generalization of the quasiparticle interaction 
in Landau's theory of Fermi liquids |1U2] . An approximate solution of the Eq flow equation for 
symmetric nuclear matter and neutron matter as a function of A is shown in Fig. [37] for two different 
Fermi momenta kp (corresponding to different densities) |141] . As expected, the in- medium SRG drives 
the Hamiltonian to a form where Hartree-Fock becomes exact in the limit A — t- 0. In contrast to the 
ladder approximation based on NN-only SRG interactions evolved in free space, the same many-body 
calculation using interactions evolved with the in-medium SRG at the two-body level gives energies 
that are approximately independent of A. This indicates that truncations based on normal-ordering 
efficiently include the dominant induced many-body interactions via the density-dependent zero-, one-, 
and two-body normal-ordered terms. 

In a similar manner, the in-medium SRG can be used as an ab initio method for finite nuclei. 
Figure 1551 illustrates this with calculations of the ground-state energy of ^He |n2j . As the flow parameter 



s increases. Fig. 38(a) shows that the Eq flow and including second-order (in F) many-body perturbation 



^■^The particlc-holc contributions to the SRG flow equation for F, given by the last sum and the second line of Eq. ([5 
are neglected in the infinite matter calculations. No such approximation is made for the '^He results presented in Fig 
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Figure 38: Ground-state energy of ^He based on the in-medium SRG starting from a free -space SRG 
potential (A = 2.0 fm~^ from the N^LO potential of Ref. [20]) [142] . (a) Eq flow (solid line) and including 
second-order many-body perturbation theory (MBPT) contributions (dashed line) as a function of flow 
parameter s. (b) Convergence behavior as a function of harmonic oscillator parameter fuj with increasing 
single-particle Cmax = max(2n-|-Z). The horizontal dotted line denotes the ground-state energy obtained 
from an exact diagonalization |85j . 



theory contributions approach each other, as was the case for the infinite matter results in Fig. [371 In 



addition, the convergence behavior with increasing harmonic-oscillator spaces in Fig. 38(b) is very 
promising. Based on these preliminary calculations, the in-medium SRG truncated at the normal- 
ordered two-body level appears to give accuracies comparable to coupled-cluster calculations truncated 
at the singles and doubles (CCSD) level. Future applications to ^^O and ^°Ca will help clarify if this 
similarity to CCSD is systematic, although such a connection was suggested by White in his work on 
similar methods in quantum chemistry |143] . Finally, we note that the in-medium SRG is a promising 
method for non-perturbative calculations of valence shell-model effective interactions and operators. 

4.4 Effective operators 

As discussed in Section 12. 5[ all operators in a low-energy effective theory will evolve under a change 
of resolution. Given an initial operator that is consistent with the initial Hamiltonian (and which will 
also have many-body components), it is important to consistently evolve the operators along with the 
Hamiltonian. When this evolution is via unitary transformations, the operator evolution is well defined 
and in principle straightforward to implement. 

For Viowfc interactions, the use of Lee-Suzuki transformations provides a convenient formalism to 
evolve consistent two-body operators [SB]- In the notation of Ref. |144[ll45j . the evolution of an opera- 
tor O from a momentum cutoff Aq to lower resolution A is given by 

0(A) = -y=L=iP + ct) O (Ao) (P + u;)^=L= , (53) 

where the operator u = QuP generates the Lee-Suzuki transformation, the projection operator P 
projects onto relative momenta A; < A, and Q = 1 — P projects onto A < k < Aq. While this approach 
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is straightforward to implement at the two-body level (and has been tested), the generalization to 
higher-body components is not simple. 

For the SRG approach, the evolution with s of any other operator O is given by the same unitary 
transformation that evolves the Hamiltonian, Os = U{s)OW (s), which means that Os evolves as 

^ = 0,] = [[T,K],0,]. (54) 

Just as with the Hamiltonian Hg, this evolution will induce many-body operators even if the initial 
operator is purely two-body. If we restrict ourselves to two-body operators in relative momentum 
space, one has 



dOs{k,k') _ 2 
ds vr 



oo 







dq {e - q^) Vs{k, q) 0,(g, k') + {k'^ - q') 0,{k, q) K(g, k') . (55) 



To evolve a particular Os simultaneously with Vs., we can simply include the discretized version of 
Eq. (!55|) as additional coupled first-order differential equations. 

An alternative, more direct, approach is to construct the unitary transformation at each s explicitly, 
and then use it to transform operators. Let \ipa{s)) be an eigenstate of Hs with eigenvalue (which 
is independent of s). Then U{s) is given by 

U{s) = Y,\Ms)){i^am- (56) 

a 

In a discretized partial- wave momentum-space basis with momenta {fcj}, we can solve for the eigenvec- 
tors oi H = Hs=Q and Hs, then construct the matrix elements of U{s) [which we denote Us{ki, kj)] by 
summing over the product of momentum-space wave functions: 

Us{ki,kj) = J2 {ki\Ms)){MO)\kj) . (57) 

a 

In practice this is an efficient way to construct the unitary transformation and subsequently to evolve 
any operator in a few-body space. For example, in the harmonic-oscillator basis used for evolving 3N 
interactions (see Section [4.21) . one can simply apply Eq. (1561) in this basis. Proof-of-principle tests for 
operators in a one-dimensional model show that this works as expected [146]. 

As an example in the two-body system, we consider the momentum distribution in the deuteron |68j . 
The momentum distribution at relative momentum q is the expectation value of Aqflq (summed over 
spin substates Ms) and is proportional to the sum of the squares of the normalized S- and D-state parts 
of the deuteron wave function, M(g)^ + ty(g)^. It is not an observable (see, for example, Ref. [147] ) but is 
useful to illustrate some points. Because the SRG proceeds via unitary transformations, by construction 
no information is lost in the evolution. But even more, we can show explicitly that by evolving to a 
low-momentum interaction, we decouple the low- and high-momentum physics in a low-energy state. 



In Fig. 39(a), we plot deuteron matrix elements of a'^a^ for the Argonne Vis [IB] and CD-Bonn 
potentials, as well as for two SRG and a smooth-cutoff Vjowfc interaction evolved from Argonne fig. We 
emphasize again that matrix elements of this operator are not measurable, so one should not ask which 
of these is the "correct" momentum distribution in the deuteron; it is a potential- and scale-dependent 
quantity. It is evident that the Vs distributions have substantial momentum components only for k below 
A, and that the V^owfc distribution is very similar to the corresponding Vs distribution. Nevertheless, if 
we use the SRG- or RG-evolved operator with these deuteron wave functions, we precisely reproduce 
the momentum distribution from the original potential at all momenta and for all A. 

This result by itself is guaranteed by construction. The more interesting issue is where the strength 
of the matrix element comes from. For example, for the bare operator and the Argonne vis potential. 
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Figure 39: (a) Deuteron momentum distribution {a^a^d = 47r[n(g)^ + w(g)^] using the Argonne fig 
CD-Bonn [19] and SRG potentials evolved from Argonne to A = 1.5 fm~^ and 2fm^^ (but not 
evolving the operator), and a smooth-cutoff Viow/c interaction with A = 2fm~^ and riexp = 2. (b) Ratio 
of the deuteron momentum distribution at various momenta q evolved from the Argonne fig potential 
via the SRG to A = 1.5 fm~^ and 2fm~^, to the corresponding initial momentum distribution, as a 
function of the maximum momentum fcmax in the deuteron wave functions in the numerator. Note that 
the unevolved Argonne fig result is simply a step function at q. For comparison, we also show the result 



for a smooth-cutoff Vfowfc interaction with A = 2fm and n, 



exp 



2. For details see Ref. 



comes entirely from deuteron wave function components 
\ it is clear from Fig. 39(a) that the deuteron does not 



the momentum distribution at g = 4fm~ 
at that momentum. But when A = 2fm 
have appreciable momentum components above 2.5 fm~^ (even though V^^ has matrix elements near the 
diagonal [7]) because of decoupling. One might imagine that the operator becomes pathological to 
compensate. In fact, it is also softened. 



In Fig. 39(b) , we take the ratio of the evolved operator evaluated with the evolved wave function 
at q to the corresponding initial quantity, but include in the numerator only momenta up to fcmax- The 
numerator is thus: 



rfk / rfk>J(k;A)t/A(k,q)t/l(q,kOV^d(k';A) 



(5J 



We observe that for all g, the ratio approaches unity for large enough fcmax, as dictated by the unitary 
transformation. For A = 1.5 fm~^ or 2fm~^, larger values of g (g > 3fm~^ or g > 4fm~^, respectively) 
give results approximately independent of g, with a smooth approach to unity by /cmax ~ 1.3 A. This is 
consistent with the operator f/A(k, q) factorizing into K\{Sl) Q\{f^ for < A and q ^ X, and thus the g 
dependence cancels in the ratio ^68j. It remains to be seen whether this factorization is a general feature 
that can be understood using operator product expansion ideas [5 ^1146] . The immediate point is that 
the flow of the operator strength weighted by ipd is toward lower momenta (softened). Furthermore, 
there is no sign of fine-tuning in the evolved operator, because even a rough variational Ansatz for the 
wave function is accurate. More work on operators is needed, but thus far these observations seem to 
be general |146] . 
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5 Applications to infinite matter 



The physics of nucleonic matter ranges over exciting extremes: from universal properties at low densi- 
ties [M1I148] that can be probed in experiments with ultracold atoms [149] : to understanding nuclear den- 
sities in known and new exotic nuclei [ETlllSUj : to dense matter in neutron stars and supernovae |151[|152j . 
At very low densities, the details of pion exchanges are not resolved and pionless EFT is well suited for 
capturing large scattering-length physics with systematic improvements from including effective range 
and higher-order terms [30l[32]. For very dilute neutron matter, there are no length scales associated 
with the interaction, because the scattering length is large fcpas ^1- As a result, only the density 
sets the scale and the properties of neutron matter are universal |33j0 With increasing density the 
effective range correction weakens the interactions of nucleons with higher momenta. This simplifies 
the many-body problem and has been used to predict the energy of neutron matter with controlled 
uncertainties for densities p ~ po/10 ISj (with saturation density po = 0.16 fm~'^). 

For nuclear densities, Pauli blocking suppresses the nonperturbative physics due to weakly and 
nearly bound states and the large scattering lengths are dissolved [9l ll54l[T55] . The difficulties in the 
theory of nuclear matter arise from the strong short-range repulsion and strong short-range tensor forces 
in phenomenological or chiral EFT potentials. The resulting coupling between low and high momenta is 
still present in the G matrix, as discussed in Section 13.51 This renders the Bethe-Brueckner-Goldstone 
expansion nonperturbative |156j . As detailed in Sects. [2111 and 12. 4( the RG evolution softens these short- 
range parts, which makes the many-body calculation more controlled. In this section, we review the 
application of low-momentum interactions to infinite matter with p ~ po and highlight the possibility 
of perturbative approaches [9| [T6l, [T7] . The findings can provide a starting point to develop a new power 
counting for nuclear matter. We focus on calculations including 3N forces and so only mention the work 
of Refs. [TFmT63] . 

5.1 Nuclear matter 

Over the last fifty years, an accurate prediction of nuclear matter starting from nuclear forces has been 
a theoretical milestone on the way to finite nuclei, but has proved to be an elusive target. Despite 
the long-term emphasis on the infinite uniform system, most advances in microscopic nuclear structure 
theory over the last decade have been through expanding the reach of few-body calculations. This 
has unambiguously established the quantitative role of 3N forces for light nuclei {A ^ 12) p2l - [2^ [57 | 
[60l[6T]. However, until recently few-body fits have not sufficiently constrained 3N force contributions 
at higher density such that nuclear matter calculations are predictive. Nuclear matter saturation is 
very delicate, with the binding energy resulting from cancellations of much larger potential and kinetic 
energy contributions. When a quantitative reproduction of empirical saturation properties has been 
obtained, it was imposed by hand through adjusting phenomenological short-range three-body forces 
(see, for example, Refs. [5311164] ). 

Progress for controlled nuclear matter calculations has long been hindered by the difficulty of the 
nuclear many-body problem when conventional nuclear potentials are used. Recent calculations over- 
come these hurdles by combining systematic starting Hamiltonians based on chiral EFT j^OtlH] with RG 
methods [5l[6] to soften the short-range repulsion and short-range tensor components of the initial chiral 
interactions [12]. The RG evolution to low momenta leads to contributions in the particle-particle chan- 
nel that are well converged at second order in the potential, suggesting that perturbative approaches can 
be used in place of the Bethe-Brueckner-Goldstone hole-line expansion [^[TB]. The Weinberg eigenvalue 
analysis discussed in Section 12.41 provides quantitative backing to these observations [T1[T2] . 

^''^In this regime, nuclear matter exists as a mixture of nucleons and nuclei that collapses to clusters at low tempera- 
tures [153) . 
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Figure 40: Nuclear matter energy per particle as a function of Fermi momentum kp at the Hartree-Fock 
level (left) and including second-order (middle) and particle-particle-ladder contributions (right), based 
on evolved N^LO NN potentials and 3N forces fit to Eau and r4He [IS]. Theoretical uncertainties are 
estimated by the NN (lines) and 3N (band) cutoff variations. 



The key nuclear matter results are summarized in Fig. HQl which shows the energy per particle of 
symmetric nuclear matter as a function of Fermi momentum kp, or the density p = 2fc|/(37r^). A 
grey square representing the empirical saturation point reflects the ranges of nuclear matter saturation 
properties predicted by phenomenological Skyrme energy functionals that most accurately reproduce 
properties of finite nuclei. Nuclear matter is calculated in three approximations: Hartree-Fock (left) 
and including approximate second-order (middle) and summing particle-particle-ladder contributions 
(right). These are the first results for nuclear matter based on chiral NN and 3N interactions. The 
technical details are given in Refs. [9l|T^ and work is in progress to improve the 3N treatment. 

The calculations of Fig. gO] start from the N^LO NN potential (A = 500 Me V) of Ref. pO]. This 
NN potential is RG-evolved to smooth-cutoff low- momentum interactions Viowfc using the techniques 
of Section 13.21 Based on the universality of Viowfc discussed in Section 12. 3[ we do not expect large 
differences starting from different N^LO potentials. The N^LO 3N forces of Section 14.11 are taken as 
a truncated basis for low-momentum 3N interactions]^ where the cd and ce couplings have been fit 
for various cutoffs to the '^H binding energy and the ^He matter radius [16]. (In the future, it will be 
possible to include consistently evolved three-body forces starting from chiral EFT using the recent 
advances in extending the SRG methods beyond the NN level [5^1137^1138] .) The 3N force fits to the 
^He radius improve the cutoff independence significantly compared to fitting to A = 3,4 energies only 
(see Fig. 6 in Ref. [9]). The calculations of Fig. HO] use the same 3N regulator of Ref. |127yi28] . but with 
a 3N cutoff A3NF that is allowed to vary independently of the NN cutoff. The shaded regions in Fig. 140) 
show the range of results for 2.0 fm~^ < A3NF < 2.5 fm~^ at fixed A = 2.0 fm~^. These predictions are 
particularly sensitive to uncertainties in the q coefficients (also for neutron matter, see Section [5.2p . 
which raises the possibility of using nuclear matter to constrain some of the q couplings. 

The Hartree-Fock results show that nuclear matter is bound even at this simplest level. A calculation 
without approximations should be independent of the cutoffs, so the spread in Fig. HQ] sets the scale 
for omitted many-body contributions. The second-order results show a dramatic narrowing of this 
spread, with predicted saturation consistent with the empirical range. The narrowing happens across 
the full density range. This is strong evidence that these encouraging results are not fortuitous, but that 
cutoff independence has been reached at the level of 1-2 MeV per particle. The controlled theoretical 

^^This assumes that the coefficients of the long-range 27r-exchange part of 3N forces are not modified by the RG 
evolution. 
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Figure 41: (a) Nuclear matter energy of Fig. HOjat the particle-particle-ladders level compared to SRG- 
evolved chiral NN interactions, (b) Nuclear matter energy of Fig. HQ] at the particle-particle-ladders 
level compared to NN-only results for two representative NN cutoffs and a fixed 3N cutoff. 



uncertainties come as a result of the combination of chiral EFT and RG methods. For all cases, the 
compressibility K = 190-240 MeV is in the empirical range. To our knowledge, these are the first 
nuclear forces fit only to A ^ 4 nuclei that predict realistic saturation properties. Similar nuclear 
matter energies are found in Fig. 41(a) when the SRG is used to evolve the NN potential. This helps 
support the general nature of the 3N force fit. On the other hand, the difference of 2 MeV per particle 
at saturation and above enlarges the theoretical uncertainty. 

The decreases in cutoff dependence in Fig. HQ] with more complete approximations is necessary but 
not sufficient to conclude that the calculations are under control. Indeed, approximations that are 
independent of the cutoff will shift the answer but not widen the error band from cutoff variation. 
The theoretical errors arise from truncations in the initial chiral EFT Hamiltonian, the many-body 
approximations, and the approximation of 3N forces. The latter is particularly uncertain here because 
it involves long-range contributions that are independent of the cutoff. Many-body corrections to 
the current approximations include higher-order terms in the hole-line expansion and particle-hole 
corrections. An approach such as coupled cluster theory that can perform a high-level resummation 
will be necessary for a robust validation. 

The evolution of the cutoff A to smaller values is accompanied by a shift of physics. In particular, 
effects due to iterated tensor interactions, which peak in the relative momentum range k ~ 4fm~^ (and 
thus lead to saturation at too high density), are replaced by three-body contributions. We emphasize 
that the relative importance of contributions to observables from the tensor force or from 3N forces 
are scale or resolution dependent. Renormalizing the potential to lower momenta redistributes these 
contributions and makes the many-body approach tractable. The role of 3N forces for saturation is 



demonstrated in Fig. 41(b) The two pairs of curves show the differences between the nuclear matter 
results for NN-only and NN plus 3N interactions. It is evident that saturation is driven by 3N forces (9] 
[16] . Even for A = 2.8 fm~^, which is similar to the lower cutoffs in chiral EFT potentials, saturation is at 
too high density without 3N forces. While 3N forces drive saturation for low-momentum interactions, the 
3N contributions are not unnaturally large. This is demonstrated in Fig. |l2]by the ratio of (Vsn) / {V\owk) 
for the triton, alpha particle, and nuclear matter at various densities. 
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Figure 42: Ratio of V3N to l^owfc expectation values. For nuclear matter, the expectation values are 
obtained by the Feynman-Hellman method and (Vsn) corresponds to the total 3N contribution at second 
order from Ref. |9], whereas the values shown for and ^He are the largest 3N contributions from 
Ref. (because of cancellations in the light nuclei). 



In Fig. HO] the particle-particle-ladder sum is little changed from second order except at the lowest 
densities shown, where the presence of two-body bound and nearly bound states necessitate a nonper- 
turbative summation. The rapid convergence for p ~ po ^lay justify in part the application of chiral 
EFT to nuclear matter in perturbation theory |165til68j . but a detailed study is needed, in particular 
regarding the cutoff dependence and the contributions from contact interactions and 3N forces beyond 
A excitations. 

While nuclear matter has lost its status to light nuclei as the first step to nuclear structure, it is still 
key as a step to heavier nuclei. The low-momentum results, while not conclusive, open the door to ab 
initio density functional theory based on expanding about nuclear matter [67 11169] . This is discussed 
further in Section 16.41 



5.2 Neutron matter 

Neutron matter provides a different perspective from symmetric nuclear matter, because only the long- 
range 27r-exchange ci and C3 terms of N^LO 3N forces contribute [17 11170] . Density-dependent two-body 
interactions V^^ based on these terms for general momentum and spin configurations and including all 
exchange diagrams were derived in Ref. p7] by summing the third particle over occupied states in the 
Fermi sea. This corresponds to the normal-ordered two-body part of 3N forces. Effective interactions of 
this type have been studied in the past using 3N potential models and approximate treatments (see, for 
example, Refs. [171(1172] ) and recently in the framework of in-medium chiral perturbation theory [173] . 
The partial- wave matrix elements in Fig. H3] show that the resulting V-^^ is dominated in neutron matter 
by a repulsive central part. The density dependence of the two-body matrix elements depends on the 
partial wave. In the ^Sq channel this dependence can be approximately parameterized by a power of the 
Fermi momentum. In the spin-triplet channels this is more complex due to the different momentum and 
density dependences of the various operator structures in Vs^. In addition to the density dependence, 
in general V^3n depends on the two-body center-of-mass momentum, but this dependence is weak and 
the P = approximation for V^sn leads to energies that agree very well with the exact result at the 
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Figure 43: Diagonal momentum-space matrix elements of the density-dependent two-body interaction 
l^3N for P = in the ^Sq and the spin-triplet P-wave channels. Results with A3NF = 2.0 fm~^ are shown 
versus relative momentum k for different Fermi momenta kp = 1.0, 1.2, 1.4 and 1.6 fm~^ (increasing in 
strength). For the P-wave channels and kp = 1.6 fm~^, the dotted lines represent the central parts 
(degenerate in J) of I^sn, whereas the dashed lines include the central plus tensor interactions. For 
details see Ref. [T7]. 

Hartree-Fock level [T7] . 

Using T^aN simplifies the neutron matter calculation based on a loop expansion around the Hartree- 
Fock energy [T7]. Figure HH shows the energy of neutron matter as a function of the density p = /c|/(37r^) 
at second-order. These are the first results for neutron matter based on chiral EFT interactions and 
including N^LO 3N forces. As discussed for nuclear matter, the RG evolution of N^LO NN potentials to 
low momenta renders the many-body calculation more controlled. The results of Ref. [T7] suggest that 
neutron matter is perturbative at nuclear densities. This is based on small second-order contributions, 
self-energy corrections being negligible, and a generally weak cutoff dependence in Fig. 44(a) For 
low-momentum interactions, N^LO 3N forces provide a repulsive contribution to the energy due to the 
repulsive central part in V^^. 



Figure 44(b) shows that the theoretical uncertainties of the neutron matter energy are dominated 
by the uncertainties in the Cj coefficients, in particular the C3 part. The Ci and C3 variation leads to an 
energy uncertainty of ±1.5 MeV per particle at saturation density. This is larger than the uncertainties 
probed by cutoff variations, which are due to the many-body approximations or from neglected short- 
range many-body interactions. We note that other recent neutron matter calculations lie within the 
Ci uncertainty band of Fig. 44(b) The spread in the energy band implies microscopic constraints for 
the symmetry energy and its density dependence [T7|. The resulting range for the symmetry energy 
04 = (31.6-34.7) MeV (where the error is dominated by the uncertainty of C3) is very useful given 
the scale of the empirical range = (25-35) MeV |174j . Other results for neutron matter with low- 
momentum interactions include studies of the equation of state at finite temperature [170] and of 
quasiparticle interactions and Fermi liquid parameters using a particle-hole RG approach [175H178j . 



5.3 Pairing 

Superfluidity plays a central role in strongly-interacting many-body systems. Pairing in infinite matter 
impacts the cooling of isolated neutron stars |179H181] and of the neutron star crust following X-ray 
bursts in accreting neutron stars |182H184j . In addition, studies of superfluidity are used to constrain 
and improve calculations of pairing in nuclei (see, for example, Ref. [97]). This section provides a brief 
review of pairing in infinite matter based on low-momentum interactions. 

shows superfiuid pairing gaps in neutron matter, obtained by solving the BCS gap 



Figure 45 (a 



equation with a free spectrum. At low densities (in the crust of neutron stars), neutrons form a ^Sq 
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Figure 44: (a) Energy per particle of neutron matter as a function of density p at second order. The 
results are based on evolved N^LO NN potentials and N^LO 3N forces. Theoretical uncertainties 
are estimated by varying the NN cutoff (lines) and the 3N cutoff (band for fixed A = 2.0 fm~^). 
(b) Theoretical uncertainties of the second-order energy due to the uncertainties in the c\ and C3 
coefficients of 3N forces. The resulting band is compared to other neutron matter results, for details 
see Ref. [H]. 



superfluid. At higher densities, the S-wave interaction is repulsive and neutrons pair in the '^P2 channel 



(with a small coupling to due to the tensor force). Figure 45(a) demonstrates that the ^Sq BCS 
gap is practically independent of nuclear interactions, and therefore strongly constrained by NN phase 
shifts [96]. This includes a very weak cutoff dependence for low-momentum interactions with sharp 
or sufficiently narrow smooth regulators with A > 1.6 fm~^. The inclusion of N^LO 3N forces leads to a 
reduction of the ^Sq BCS gap for Fermi momenta k-p > 0.6 fm~ [17J. This reduction becomes significant 
for densities where the gap is decreasing and agrees qualitatively with results based on 3N potential 
models (see, for example, Ref. |185j ). Two-nucleon interactions are well constrained by scattering data 
for relative momenta k < 2fm~^ [5]. The model dependencies at higher momenta show up prominently 
in Fig. 45(a) in the ^P2— ^F2 gaps for Fermi momenta kp > 2fm~^ |186j . 



For pairing in nuclei, calculations with V^owfc as pairing interaction provide a good starting description 
of neutron- neutron and proton-proton gaps [7 ^[501[^I188] . as discussed in Section This is very 
promising and also consistent with comparisons of the ^Sq BCS gaps in neutron matter with results based 
on the empirical Gogny interaction |189[I190] . In nuclei, however, partial waves beyond the standard 
^So channel can contribute to the pairing interaction and play an interesting role for the T = 1, J = 
pair formation (because paired nucleons are not in back-to-back-momentum configurations) |191j . 

Understanding many-body effects beyond the BCS level constitutes an important open problem. 
For recent progress and a survey of results, see for instance Ref. |192] . Here we review approaches to 
this problem using the RG. At low densities, induced interactions due to particle-hole screening and 
vertex corrections are significant even in the perturbative kpa limit |193[I194] and lead to a reduction 
of the S-wave gap by a factor (4e)~^/^ ~ 0.45, 
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Figure 45: (a) The ^Sq (left panel) and ^P2— '^F2 (right panel) superfluid pairing gaps A at the Fermi 
surface as a function of Fermi momentum kp in neutron matter. The gaps are obtained from charge- 



dependent NN interactions at the BCS level. For details see Refs. [9611186] . (b) Top panel: Comparison 
of the ^So BCS gap to the results including polarization effects through the phRG, for details see 
Ref. [176j . and to the results of Wambach et al. [IHT]. Lower panel: Comparison of the full superfluid 
gap A to the BCS gap Aq and to the Fermi energy e-p. 



This reduction is due to spin fluctuations, which are repulsive for spin-singlet pairing and overwhelm 
attractive density fiuctuations0 

Following Shankar |197j . Ref. |176j developed a nonperturbative RG approach for neutron mat- 
ter, where induced interactions are generated by integrating out modes away from the Fermi surface. 
Starting from l^owfc? the solution to the RG equations in the particle-hole channels ("phRG") includes 
contributions from successive particle-hole momentum shells. The phRG builds up many-body effects 
similar to the two-body parquet equations, and efficiently includes induced interactions to low-lying 
states in the vicinity of the Fermi surface beyond a perturbative calculation. The phRG results |176] 
for the ^So gap are shown in Fig. 45(b), where induced interactions lead to a factor 3-4 reduction to 
a maximal gap A 0.8 MeV. For the lower densities, the phRG is consistent with the dilute result0 
A/Ao = (4e)~^/^, and at the larger densities the dotted band indicates the uncertainty due to an 
approximate self-energy treatment. 

Non-central spin-orbit and tensor interactions are crucial for ^P2— ^F2 superfluidity. Without a 
spin-orbit interaction, neutrons would form a ^Pq superfluid instead. The flrst perturbative calculation 
of non-central induced interactions shows that ^P2 gaps below 10 keV are possible (while second-order 
contributions to the pairing interaction are not substantial |^nd/^owfc| < 0.5) |177] . This arises from a 
repulsive induced spin-orbit interaction due to the mixing with the stronger spin-spin interaction. As 
a result, neutron P-wave superfluidity (in the interior of neutron stars) may be reduced considerably 
below earlier estimates. This implies that low-mass neutron stars cool slowly [179H181] . Smaller values 
for the ^P2 gap compared to Fig. 45(b) are also required for consistency with observations in a minimal 
cooling scenario [198j . 



^^In finite systems, the spin and density response differs. In nuclei witfi cores, tire low-lying response is due to surface 
vibrations. Consequently, induced interactions may be attractive, because the spin response is weaker |1951ll96j . 

"'^''For fcp ~ 0.4 fm~^, neutron matter is close to the universal regime, but theoretically simpler due to an appreciable 
effective range k-pr^ « 1 [33]. 
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Figure 46: Ground-state energy of ®Li versus oscillator parameter HQ for different SRG-evolved inter- 
actions with A = 3.0,2.0,1.5 and 1.0 fm~^. The initial interaction is the N'^LO NN-only potential of 
Ref. [20]. The NCSM results clearly show improved convergence with the maximum number of oscil- 
lator quanta A^max for lower cutoffs. Since 3N interactions are neglected, the different NN calculations 
converge to different ground-state energies. For details see Ref. [H]. 

6 Applications to finite nuclei 

The advantages of low-momentum interactions discussed in Section [5] for infinite matter calculations 
also apply to finite nuclei. At lower resolution, perturbative methods generally become more effective, 
finite-basis expansions converge more rapidly, and short-range correlations are substantially weakened. 
Furthermore, the freedom to vary the cutoff with the RG provides a useful tool to diagnose truncation 
errors and assess the quality of many-body approximations. In this section, we review applications of 
low-momentum interactions to finite nuclei, including ab initio calculations of light and medium-mass 
nuclei, shell model calculations, scattering and reactions, and energy- density functional approaches. 

6.1 Ab initio calculations 

The no-core shell model (NCSM) is a versatile method to calculate ground-state energies and spectra 
of light nuclei [65]. In Ref. [H], the NCSM was used to study the convergence of nuclei up to ^Li using 
RG-evolved chiral NN interactions. Because no additional Lee-Suzuki transformations are required, the 
energy calculations satisfy the variational principle for a given low-momentum Hamiltonian. As shown 
in Fig. HHl the convergence of ground-state energies improves rapidly for SRG potentials with decreasing 
A, with similar convergence improvements for Viowfc interactions. 

When the RG equations are truncated at the two-body level, the evolution is only approximately 
unitary and converged energies for A ^ 3 vary with A or A. This approximation is systematic, however, 
and for useful cutoff ranges, the energy variation shown in Fig. HT] is comparable to natural-size trun- 
cation errors in chiral EFT, with no unnaturally large contributions from omitted three-body forces for 
these light nuclei. Further evidence that the evolution preserves the hierarchy of the underlying chiral 
EFT is shown in Figs. |35] and |36l 

Extensions of ab initio methods to heavier and neutron-rich nuclei are a frontier of nuclear theory. 
Coupled-cluster (CC) theory is the prime method for systems with up to 100 electrons in quantum 
chemistry |200j and a powerful method for nuclei for which a closed-shell reference state provides a 
good starting point [86 tl201] . For '^H and ^He, CC results agree with the corresponding Faddeev and 
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Figure 47: Ground-state energies of (a) ^He, (b) ^Li, and (c) ^Li as a function of A for SRG interactions 
evolved from the N^LO NN-only potential of Ref. [20]. Error bars for larger A values are from extrapo- 
lations in A^max- The arrow marks the experimental energy. The characteristic increase in Egg at small 
A signals the modification of the long-range attractive interaction [138] . For details see Ref. 
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Figure 48: CCSD ground-state energy of ^°Ca as function of the oscillator frequency with increasing 
basis size N, based on (a) an RG-evolved low-momentum NN interaction [13j and (b) a N^LO NN 
potential [199] . The difference highlights the importance of 3N interactions for ground-state energies. 



Faddeev-Yakubovsky energies [_13j. Combined with rapid convergence for low-momentum interactions, 
CC theory has pushed the limits of accurate calculations to medium-mass nuclei and set new benchmarks 
for ^^O and ^°Ca [T3|. Using an angular-momentum-coupled scheme, it is possible to extend CC theory 
to very large spaces (15 major shells on a single processor) and to obtain near-converged ground-state 
energies for spherical nuclei, ^°Ca, ^*Ca, and ^^Ni, based on a N^LO NN potential |199j . The CC 
developments for medium-mass nuclei are shown in Fig. HHl where the critical importance of 3N forces 
for ground-state energies is evident. 

Recently, a combination of nuclear and atomic physics techniques led to the first precision mea- 
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Figure 49: Convergence as a function of the basis size for the HH ^He and ^He ground-state energies 
and for the CCSD ground-state energy of ^He based on chiral low-momentum NN interactions Viowfc for 
a range of cutoffs A = 1.8, 2.0 and 2.4 fm'^ For details see Ref. [15]. 
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Figure 50: Relative contributions \AE/E\ to the binding energy of ^He at the CCSD level from Viowfc 
as well as normal-ordered zero-, one-, two-body and residual three-body parts of 3N forces [202] . 



surements of masses and charge radii of the helium halo nuclei, ®He |203[l204j and ^He [20511206] . with 
two or four weakly-bound neutrons forming an extended halo around the ^He core. In Fig. HHl results 
are shown for the ground-state energies of helium nuclei based on chiral low-momentum NN interac- 
tions [15]. This combines the RG evolution with the exact hyperspherical-harmonics (HH) expansions 
for ^He and CC theory for '^He (see also Ref. [86]), which have the correct asymptotic behavior of the 
wave function. The cutoff variation in Fig. |19] highlights the importance of 3N interactions. For all 
studied cutoffs, the NN-only results underbind ^He [T5l|86]. Therefore, the helium isotopes probe 3N 
effects beyond the overall repulsion observed in infinite nuclear and neutron matter. 

Finally, in Fig. [50] the first CC results with 3N forces show that low-momentum 3N interactions are 
accurately treated as normal-ordered zero-, one-, and two-body terms, and that residual 3N forces can 
be neglected [202] . This is very promising for developing tractable approximations to handle many-body 
interactions (see Section 1^731) . and supports the idea that phenomenological adjustments in shell model 
interactions ( "monopole shifts" ) are due to 3N contributions [207] . These monopole shifts are enhanced 
in neutron-rich nuclei and have a pivotal impact on shell closures and the location of the neutron drip 
line, as discussed in Section 16^31 
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Figure 51: (a) Neutron analyzing power Ay for neutron-deuteron scattering at £"„ = 3 MeV and 
^cm = 104° as a function of RG cutoff A (left panel) and SRG parameter A (right panel). The horizontal 
line at 10^ Ay = 5.86 is the experimental value from Ref. |208] . (b) Maximum of the neutron analyzing 
power Ay for neutron-^H scattering at En = 3.5 MeV as function of RG cutoff A (left panel) and SRG 
parameter A (right panel). For details see Ref. |209] . 



6.2 Scattering and reactions 



In order to probe the sensitivity of 3N and 4N scattering observables to short-range physics, Ref. |209 



used low-momentum NN interactions evolved with the RG or SRG to study the variation with the cutoff 
A or A. The cutoff dependence due to truncating the RG equations at the two-body level provides a 
measure of the sensitivity to neglected short-range 3N (and higher-body) forces. Comparing the results 
shown in Fig. [5T]one notices that the cutoff dependence of 3N scattering observables is much weaker than 



the one observed for 4N observables. This weak cutoff dependence in Fig. 51(a) implies that short-range 
3N forces are not likely to solve the long-standing discrepancies of the nucleon-deuteron Ay in elastic 
scattering with data. This is indeed what has been found when short-range or the leading long-range 3N 
interactions are included. Nucleon-deuteron Ay in elastic scattering and the space star differential cross 
section for breakup barely change by adding a 27r-exchange 3N force |210H212] . an effective 3N force 
due to the explicit A-isobar excitation |213H216] . or N^LO 3N forces [128] . In contrast, the maximum 
of the neutron analyzing power Ay for neutron-^H scattering seems to be more sensitive to short-range 



many-body forces as demonstrated by the more pronounced dependence on the cutoff in Fig. 51(b) 
This is consistent with previous findings |217] obtained with various NN potentials. 

Neutrinos play a crucial role for the physics of stellar collapse, supernova explosions and neu- 
tron stars |218U219j . The first calculations |220[l221j of neutrino processes based on chiral EFT 
and RG-evolved interactions have focused on neutrino reactions involving two nucleons: neutrino-pair 
bremsstrahlung and absorption, NN -H- NNuV, which are key for muon and tau neutrino production 
in supernovae, and neutrino inelastic scattering, vNN vNN . These processes are determined by 
the spin relaxation rate I/to- using a unified approach to neutrino interactions in nucleon matter 



In supernova simulations, the standard rates for bremsstrahlung are based on the one-pion exchange 
approximation to NN interactions |222] . Calculations based on chiral EFT and RG-evolved interactions 
show that shorter-range non-central forces reduce the spin relaxation rate and therefore the rate for 
bremsstrahlung significantly |22Utl221] . In addition, the spin relaxation rate sets the scale for energy 
transfer in inelastic scattering, which is presently not included in simulations. The spin relaxation rates 
of Ref. |221] lead to a mean-square energy transfer that is comparable to the incoming neutrino energy, 
indicating that NN collisions may be important for energy transfer in supernovae. 

Quaglioni and Navratil have recently developed an ab initio approach to nuclear reactions based 
on merging the No-Core Shell Model (NCSM) with the Resonating Group Method (RGM), where the 
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Figure 52: (a) No-Core Shell Model (NCSM) and Resonating Group Method (RGM) applied with an 
SRG-evolved N^LO potential (NN-only, A = 2.02 fm~^) to calculate nucleon-alpha elastic scattering 
phase shifts at En = 17MeV |223j . (b) Comparisons are made to the differential cross section (bottom 
panel) and the analyzing power (top panel) from Karlsruhe polarized-neutron experiments [223]. 

many-body Hilbert space is spanned by cluster wave functions describing a system of two or more 
clusters in relative motion [22411225] . Figure [52] shows nucleon-alpha elastic scattering at En = 17MeV 
using an SRG-evolved N^LO potential [223 J. As with the ground-state NCSM calculations, the use 
of low- momentum interactions in the NCSM/RGM description of nucleon-nucleus scattering results in 
faster convergence with respect to the basis size. This is a very promising start. 

6.3 Shell-model approaches 

Shell-model applications were one of the main motivations for developing low-momentum interac- 
tions [226(1227] . The practical benefit of low- momentum interactions in shell model applications is 
shown in Fig |53l where we compare S-wave relative harmonic-oscillator matrix elements for a low- 
momentum cutoff A = 2.0 fm~^ and a 1.0 GeV cutoff A = 5.0 fm^^. We find that for lower cutoffs 
the matrix elements decrease quite rapidly and become small for \n — n'\ ~ 5 — 10. This is not the 
case for interactions with high-momentum components, which require basis states up to ~ 50 shells for 
convergence. 

In conventional approaches to shell-model effective interactions, the strong short-range repulsion is 
tamed by performing a ladder resummation of the NN potential to obtain a G matrix [HI]. However, 
the G matrix resummation introduces an inconvenient starting-energy dependence and requires further 
approximations in practice. Furthermore, as discussed in Section 13. 5[ G matrices do not decouple 
low- and high-energy states. As shown in Fig. [SU low- lying spectra of ^^O and ^^^Te based on low- 
momentum NN interactions are of comparable quality to using G matrix resummations that are difficult 
to generalize to 3N forces and involve starting-energy dependences. In addition, as a consequence of 
decoupling, perturbative valence shell-model calculations starting from low-momentum interactions are 
under better control than the corresponding G matrix calculations. We refer to Ref. [66j for a review 
on low-momentum interactions in perturbative calculations of valence shell-model effective interactions. 
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Figure 53: Relative harmonic-oscillator matrix elements {n' = Ol'\Vi^^j^\nl) versus radial quantum 
number n for a low-momentum cutoff A = 2.0 fm~^ and a 1.0 GeV cutoff A = 5.0 fm^^. Results are 
shown for S-wave matrix elements. In both cases, Viowfc is evolved from the Argonne vis potential and 
the harmonic oscillator parameter is /lo; = 14 MeV |132] . 
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Figure 54: Low-lying spectra using shell- model effective Hamiltonians derived from Viowfc and G matrix 
interactions starting from the CD-Bonn potential [19] for (a) ^'^O and (b) ^^^Te [66] . 



Moreover, there are promising applications of low-momentum interactions in the Gamow shell model 
to handle continuum states [228f[230j . 

The neutron drip-line, which is the limit of neutron-rich nuclei, evolves regularly from light to 
medium-mass nuclei except for a striking anomaly in the oxygen isotopes. This anomaly is not repro- 
duced in theories derived from two-nucleon forces. In Ref. [231] . the first microscopic explanation of the 
oxygen anomaly based on low-momentum 3N forces was presented. As shown in Fig. [551 the inclusion 
of 3N interactions at N^LO or due to A excitations leads to repulsive contributions to the interactions 
among valence neutrons that change the location of the neutron drip-line from ^^O to the experimentally 
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Figure 55: Ground-state energies of neutron-rich oxygen isotopes measured from the energy of ^^O. The 
experimental energies of the bound oxygen isotopes i6-24q ^^^^ included for comparison. The left panel 
(a) shows the energies obtained from the phenomenological forces SDPF-M |232y233j and USD-B |234] . 
The middle panel (b) gives the energies obtained from a G matrix and including Fujita-Miyazawa 3N 
forces due to A excitations [231] . The right panel (c) presents the energies calculated from V^owfc and 
including chiral EFT 3N interactions at N^LO as well as only due to A excitations [231] . The changes 
due to 3N forces based on A excitations are highhghted by the shaded areas. 



observed ^^O. This 3N mechanism is robust and general, and therefore expected to impact predictions 
of the most neutron-rich nuclei and the synthesis of heavy elements in neutron-rich environments. 



6.4 Density functional theory 

While most advances in microscopic nuclear structure theory over the last decade have been through 
expanding the reach of few-body calculations, infinite nuclear matter is still a key step to heavier nuclei. 
In particular, the promising results using low-momentum interactions open the door to ab initio density 
functional theory (DFT) both directly (through orbital-based methods) and based on expanding about 
nuclear matter |169j . This is analogous to the application of DFT in quantum chemistry and condensed 
matter starting with the uniform electron gas in local-density approximations and adding constrained 
derivative corrections. Phenomenological energy functionals (such as Skyrme) for nuclei have impressive 
successes but lack a (quantitative) microscopic foundation based on nuclear forces and seem to have 
reached the limits of improvement with the current form of functionals |235ll236j . Furthermore, without 
theoretical understanding of errors, extrapolations to the limits of nuclear binding are uncontrolled. 

A recent review [HZj describes various theoretical paths to ab initio (meaning microscopic) DFT. One 
formal constructive framework for Kohn-Sham DFT is based on effective actions of composite operators 
using the inversion method |237H242] . This is an organization of the many-body problem that is based 
on calculating the response of a finite system to external, static sources rather than seeking the many- 
body wave function. Other paths include many-body perturbation theory or RG methods applied to 
effective actions p7 tl243[r244] . Either way necessitates a tractable expansion that is controllable in the 
presence of inhomogeneous sources or single-particle potentials. This is problematic for conventional 
nuclear forces, for which the single-particle potential needs to be tuned to enhance the convergence of 
the hole- line expansion |123[ll56j . but is ideally suited for low- momentum interactions. 

The nuclear matter results of Fig. HD] and Ref . [HlIIS] imply that exchange correlations are tractable, 
and this motivates the derivation of an ab initio density functional based on low-momentum interactions, 
for example, using a density matrix expansion |169j (see also Ref. [33 |I245] ). This is one of the goals 
of the SciDAC universal nuclear energy density functional project (UNEDF) [18]. Future work can use 
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Figure 56: Neutron/proton pairing gaps (upper/lower panel) obtained from calculated and experimental 
odd-even mass differences along isotopic/isotonic chains. The calculated gaps use the ^Sq part of V^owA; 
as pairing interaction in density functional calculations (plus Coulomb for proton gaps) [188(1246] . 

EFT and RG to identify new terms in the functional, to quantify theoretical errors for extrapolations, 
and to benchmark with ab initio methods (such as with coupled-cluster theory for medium-mass nuclei). 

The theoretical errors of nuclear matter calculations based on low-momentum interactions, while 
impressively small on the scale of the potential energy per particle, are far too large to be quantitatively 
competitive with existing energy functionals. However, there is the possibility of fine tuning to heavy 
nuclei, of using EFT/RG to guide next-generation functional forms, and of benchmarking with ab 
initio methods for low-momentum interactions. Work in these directions is in progress (see for example 
Ref. [SZ])- Overall, these results are very promising for a unified description of all nuclei and nuclear 
matter but much work is left to be done. 

Recent results using existing energy density functionals in the particle-hole channel (to build a 
reasonable self-consistent single-particle basis) combined with low- momentum interactions Vio^k in the 
pairing channel suggest that an ab initio DFT treatment of pairing is feasible [75 | [5I HIWlll88lll91j . 
Figure [56] shows that the resulting neutron-neutron and proton-proton pairing gaps in semi- magic nuclei 
are remarkably consistent experiment [188[I246] . Systematic investigations of theoretical corrections, 
including 3N forces, are underway. 

7 Summary and outlook 

As illustrated in the preceding sections, the use of low-momentum interactions has positively impacted 
few- and many-body calculations for nuclear structure and reactions in all regions of the nuclear chart. 
The advantages of softened potentials are unequivocal and recent work on evolving 3N forces addresses 
some of the most pressing concerns about the RG approach. There remain, however, some misunder- 
standings and misconceptions about the application and implications for nuclear problems of the RG 
and related methods and how they connect to EFT. In addition, there is much work in progress with 
associated open questions and as yet unrealized opportunities. Here we reiterate many of our main 
points by clarifying some of the misconceptions and surveying the on-going research. 



61 



— I — \ — r 
0.06 - AVI 8 with n 

0.05 - ^ ^ 



exp 



0.04 - 

■ 
■ 

D ■ 

0.03 - 

B 

0.02 - ■ 

@ d ^ ^ 1^ ^ga^ m^MMmMtma ® © @ 

^" • Binding energy [MeV] 
0-01 - ■ n D-state probability 

, ■ ♦ Asymptotic D/S ratio 

Q I T \ III \ \ \ I \ I I I I 



0.03 



0.025 



11h 



0.02 



-2.22 

-2.225 

-2.23 



0.5 



1 



2 3 4 5 

It 



10 



A [fm 1 

fa) 



0.06 - N^'lO [16] with n^^p = 6 



D 



0.05 
0.04- 
0.03- 
0.02- 
0.01- 




# # ^ ^ ^ 



# # ® © # g®€ 



# ® ® ®_ _2 225 



^ ■ • Binding energy [MeV] 
■ D-state probability 
♦ Asymptotic D/S ratio 



0.03 



0.025 



0.02 



- -2.22 



-2.23 



0.5 



1 2 
A [fm~^] 
(b) 



4 5 



Figure 57: D-state probability Pq (left axis), binding energy (lower right axis), and asymptotic 
D/S-state ratio rjd (upper right axis) of the deuteron as a function of the cutoff [6], starting from (a) the 
Argonne vis [IH] and (b) the N^LO NN potential of Ref. [20] using different smooth Viowk regulators. 
Similar results are found with SRG evolution. 



7.1 Misconceptions and clarifications 

In the following we address some of the misleading statements or common misconceptions associated 
with low-momentum interactions that we have encountered in the literature as well as informally. 

• The nuclear potential is an observable. Even if not stated explicitly, this is possibly an implicit 
belief inherited from experience (in both classical and quantum contexts) with the Coulomb po- 
tential, and the idea that the potential measures the energy of two particles at fixed separation. 
This intuition breaks down, of course, for finite-mass composite particles unless at large separa- 
tions. Because the short-range part of the potential can be modified with a unitary transformation 
without altering experimental predictions, it is not an observable by any meaningful definition. 

• Both two-nucleon phase shifts and lattice QCD calculations show that the NN interaction has a 
strongly repulsive core. We reiterate that a short-range quantum mechanical potential is not a 
measurable quantity. Fitting the change in sign of S-wave phase shifts from attractive to repulsive 
with a local potential leads to a repulsive core, but there are an infinity of phase-equivalent poten- 
tials without this feature. Similarly, the non-local interaction obtained from lattice simulations 
and converted to phase-equivalent local form for comparison to phenomenology [28] is based on a 
non-unique choice of nucleon interpolating field |247[l248j . 

• Low-momentum interactions change observable properties such as the D-state probability. As 
long as the corresponding operators are evolved together with the Hamiltonian, observables will 
remain unchanged. While it is true that the D-state probability of the deuteron is reduced 
as interactions are evolved toward lower cutoffs, this is not an observable [24911250] . Figure 1571 
illustrates some observable and non-observable properties of the deuteron. Up to numerical errors, 
the measurable binding energy and the asymptotic D/S-state ratio are independent of the low- 
momentum resolution, while the unobservable D-state probability changes dramatically as the 
short-range part of the tensor force is integrated out. 
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• Details of strong-interaction dynamics at high energy are important for low-energy nuclear struc- 
ture and reactions. More precisely, it has been claimed that a good (or at least reasonable) 
description of phase shifts above some energy (such as i^iab = 350 MeV) is required for a Hamil- 
tonian to be useful for nuclear structure. While it would seem natural that low- and high-energy 
physics are decoupled, this does not happen automatically and the misconception is often driven 
by experience with conventional NN potentials with strong short-range repulsion. The repulsion 
causes even bound states with very low energies (such as the deuteron) to have important con- 
tributions to the binding energy and other observables from high-momentum components (well 
above 2fm~^, see Fig. [T2l) of the deuteron wave function]^ For example, Ref. |251] presents cross 
sections for electron scattering from the deuteron calculated using as input a spectral function 
that is the momentum distribution times a delta function in energy. After noting the effect of 
excluding high momenta on the cross section, the conclusion is that "the data confirm the ex- 
istence of high-momentum components in the deuteron wave function". Beyond the fact that 
the short-range parts of wave functions are not observables, these conclusions reinforce the intu- 
ition that there is information in quantitatively reproducing high-energy phase shifts that is lost 
when evolving to low-momentum interactions. However, due to decoupling for the RG-evolved 
interactions, it is completely irrelevant if phase shifts are zero above some energy. 



Short-range correlations are important for low-energy properties of nuclei. Various electron scat- 
tering and proton scattering experiments purport to show the existence of short-range corre- 
lations in nuclei. These results are then interpreted as implying that such correlations are an 
important feature of the structure of medium to heavy nuclei (see, for example, Refs. [252(1253] ). 
As shown for the deuteron in Figs. 3(b) and 16(a), the short-range part of a wave function is 



not an observable and is radically changed by short-range unitary transformations that preserve 
all observables. Therefore any statement about short-range correlations necessarily relies on the 
(mo del- dependent) short-distance details of the Hamiltonian, which in turn depend on the reso- 
lution. As has been illustrated in this review, eliminating these correlations by evolving to low 
resolution greatly enhances the convergence of few- and many-body calculations. 



• The low-momentum RG technology fails for observables that are sensitive to high-momentum com- 
ponents. In principle, calculated observables |254[I255] from experiments such as (e, e'pN) are 
unchanged after RG evolution whether or not they involve high momentum, because the asso- 
ciated operators (which are generally not simple at high momentum) will transform to ensure 
this invariance. In practice there are many open questions, such as whether consistent operator 
evolution is difficult to approximate accurately or if factorization approximations might be more 
valid with low-momentum interactions. 



• Low-momentum interactions are not systematic. If one starts from a chiral EFT interaction 
determined to a given order in the power counting, the truncation error of the initial chiral EFT 
is exactly preserved when many-body interactions are evolved as well. If many-body interactions 
are truncated, the error remains of the same natural size for cutoffs used in practical calculations. 
Therefore, low-momentum interactions preserve the systematic nature of the initial EFT and 
EFT improvements are immediately folded into the RG evolution. In addition, the control of the 
EFT truncation errors is not lost when two-body forces are evolved exactly but the three-body 
forces are fit using chiral 3N interactions as a truncated basis. On the other hand, if one starts 
from phenomenological potentials, then it is correct that the low-momentum Hamiltonian is not 
systematically improvable because the starting point is not either. 

^®We emphasize that the need for high-momentum components in these particular calculations does not imply that the 
high-energy description is correct. 
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• Low-momentum interactions are only useful as long as calculations with "bare" chiral EFT in- 
teractions are not possible for all nuclei. We reiterate that there is nothing superior about a 
"bare" interaction as long as the RG-evolved interactions reproduce observables with a compara- 
ble truncation error. Further, converged calculations with different cutoffs are very desirable to 
quantify theoretical uncertainties. Finally, the certainty that three-body (and possibly four-body) 
forces will be needed for precision calculations implies that softer interactions will continue to be 
necessary to push to larger nuclei and to reactions even with the anticipated gains in computing 
power. (Note that normal-ordering simplifications for many-body forces are more effective for 
low- momentum Hamiltonians.) 

• There is only one Viowfc or SRG potential. In fact there are infinitely many different low-momentum 
interactions, corresponding to different initial potentials (often called "bare" interactions) and 
different choices for the evolution (cutoff value and renormalization scheme, such as sharp, smooth 
and SRG). The same holds for chiral EFT interactions, which depend on the regulator and cutoff 
values, as well as on the level of truncation. There is an analogous misconception that there is 
only one N^LO NN potential, based perhaps on the dominant use of the A = 500 MeV potential 
of Ref. [20]. Note that the full advantages of low-momentum potentials are only realized if more 
than one cutoff is used. 

• The RG evolution of nuclear forces to lower momentum improves the physics. To our knowledge, 
this is not a statement that has actually been made by a Viowfc or SRG proponent, but one which 
is often thought to be implied. However, if the initial Hamiltonian provides a poor representation 
of nature at low energies, it will not be improved by the RG evolution described here. The correct 
statement is that an exact RG evolution to lower momentum preserves the observable physics, 
while making it (generally) easier to calculate. If we could easily calculate all many-body systems 
of interest exactly with an initial Hamiltonian, then there would be no need for the RG evolution. 
(Unless changing the resolution was useful for other reasons, such as enhancing physical insight.) 

• Evolved low-momentum interactions are strongly non-local. This is misleading because it omits 
that the non-locality generated by RG evolution to lower resolution increases with the momentum 
transfer qE The local long-range (low q) part of low-momentum interactions (given by one-pion 
exchange) is therefore not modified by the evolution. Given the composite nature of nucleons, 
non-locality at short range (high q) is natural. 

• Low-momentum interactions should reproduce the initial NN-only results for A > 2. The RG 
evolution to lower resolution shifts interaction strength from two-body to few-body contributions, 
so it is guaranteed that NN-only calculations will disagree for A > 2. A related misconception is: 
Low-momentum interactions are approximations to "bare" potentials. Because induced few-body 
contributions will be truncated at some level, there will be differences from the "bare" result. 
But any "bare" interaction is also a truncated low-energy effective theory, so there should be no 
preference over low-momentum interactions with different but still natural few-body contributions. 

• Soft interactions are not realistic Hamiltonians because they overbind nuclei compared to exper- 
iment. This statement tends to follow observations that nuclei calculated with NN-only low- 
momentum interactions are increasingly overbound with increasing A. But this is only true if 3N 
forces (or the uncertainties from neglecting 3N forces) are omitted, which is not justified. The 
test in the end is whether soft Hamiltonians including many-body forces reproduce observables. 

^^This contrasts with a Pcrcy-Buck-type factorized non-local potential |256| such as V{r,r') — Viocai[(r + 
r72)](7ri/2^)-3 exp[(r - r')V/3^], for which the non-locality is independent of the momentum transfer. 
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• Many-body forces will explode after integrating out strong short-range repulsion. The calculations 
reviewed in the preceding sections demonstrate explicitly that this is not the case (for example, 
low-momentum 3N interactions are perturbative in light nuclei). Many-body contributions are 
natural in size and the decreasing hierarchy is preserved. A related misconception is: To include 
3N forces with low-momentum interactions, it is important to distinguish between "real" (or "gen- 
uine") and "induced" three-body forces. The initial and evolved 3N forces are both "effective" 
interactions. Because the EFT truncation error is preserved in the evolution, there is no physics 
reason to prefer one over the other. Each Hamiltonian has different associated many-body inter- 
actions and operators, although the long-range parts are the same as they are not modified by 
the RG evolution. 

• Weren't soft potentials ruled out for nuclear physics long ago? Soft potentials were explored 
in the mid-sixties and early seventies to improve convergence in nuclear matter and nuclei. In 
particular, separable potentials such as the Tabakin potential were considered, as well as unitary 
transformations of local, repulsive-core potentials [171159]. These NN-only soft potentials were 
abandoned primarily because they did not saturate nuclear matter in the empirical range, but at 
too-high densities and with overbinding. In contrast, with low-momentum interactions saturation 
is driven by the (natural-sized) 3N forces fit to three- and four-body properties. 

• Nuclear matter saturation with low-momentum interactions is largely due to the (long-range) two- 
pion- exchange part of 3N forces. But for conventional NN potentials, two-pion- exchange 3N forces 
are attractive for all nuclei and nuclear matter \5^\134^ . The distinction is because decoupling 
in low-momentum interactions means that 3N contributions are from the low-momentum (long- 
range) parts of 3N forces, whereas in calculations with conventional potentials, two-pion-exchange 
3N forces are also sensitive to the (model-dependent) short-range parts of wave functions. This 
leads to a net attraction for hard potentials |172j . Due to the importance of 3N forces for nuclei, 
a quantitative understanding of this change is an important open problem. 

• ^ow k is just like a G matrix. The G matrix approach leads to a softened version of an initial NN 
potential and one can see strong similarities between a G matrix and V^owfc for matrix elements 
connecting low- energy states (explicit comparisons are shown in Figs. [33] and 153]) . But the G 
matrix is an in-medium interaction that depends on the nucleus or the density and on a "starting 
energy" , while low-momentum interactions are (generally) evolved in free space. Most important 
is that the decoupling between low and high momenta that is the hallmark of low-momentum 
interactions is not achieved by the G matrix construction. The off-diagonal contributions, evident 
in Fig. |33l make V^iowfc qualitatively different from a G matrix. A consequence is that while the G 
matrix is a softened effective interaction, the off-diagonal parts will be significant in higher-order 
diagrams. This is why a perturbative expansion in the G matrix fails and something like the 
hole-line expansion is needed (with a similar G matrix resummation at the three-body level). 

7.2 Open problems and opportunities 

We finish with a survey of work in progress to extend the results described in this review together with 
open problems that may lead to new opportunities. This is, of course, not a comprehensive list, as there 
will be other efforts of which we are unaware and because new ideas constantly arise. But we hope this 
will give a reasonable vision of what to expect in the next few years. 

Chiral EFT power counting of the hierarchy of many-body forces implies that 3N forces must 
be included in almost any precision calculation of nuclear properties; this is evident already from 
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microscopic calculations of light nuclei@ Ongoing work includes the effects of 3N forces in one of three 
ways: matching to a truncated chiral EFT basis (Section I4.ip . explicit evolution with RG flow equations 
(Section I4.2p . or through normal-ordering and the in-medium SRG (Section 14. 3p . 

The recent evolution of 3N forces with SRG flow equations in a harmonic-oscillator basis will allow 
important tests of the matching approach and whether the low-momentum universality observed at the 
NN level (the collapse of different initial potentials) extends to three-body forces. These studies can 
also provide useful insights to more general EFT questions regarding the consistency of fits to NN and 
3N interactions and the power counting of three-body forces. In addition, the SRG can be extended to 
evolve 3N forces in momentum space (with the kinetic energy generator) and to generate consistent 3N 
interactions for Viowk using the SRG block diagonalization. We note that the formal equivalence between 
the block- diagonal SRG and V\owk is not yet established. Likewise, the very similar results in few- and 
many-body calculations using smooth cutoff Viowk and band-diagonal SRG (NN-only) interactions is 
not yet understood. 

The availability of low-momentum 3N interactions opens the door to many new results in ab initio 
calculations of finite nuclei using NCSM, hyperspherical harmonics, and coupled-cluster methods. An 
important open problem is to determine the range of resolutions for which the RG evolution can be used 
to soften nuclear Hamiltonians while preserving the hierarchy of many-body interactions. The study 
of residual cutoff dependencies in calculations of finite nuclei and infinite matter will help answer this 
question and provide critical tests of the RG approach. Work on the consistent evolution of operators 
will lead to new opportunities for exploring electroweak processes in nuclei and for assessing theoretical 
uncertainties of key nuclear matrix elements. 

In addition, an important direction for future work is to develop and test tractable approximations 
to handle many-body forces. The first applications of low-momentum 3N interactions to the shell model 
are very exciting and open many opportunities for exotic nuclei. The development of the in-medium 
SRG is promising to include many-body interactions through normal-ordering and to provide a non- 
perturbative approach to valence shell-model effective interactions and operators. The applications to 
reactions are another frontier. The NCSM/RGM combination has shown promise with low- momentum 
NN interactions (both Viowfc and SRG) and the extension to include 3N forces is an additional challenge. 
Revisiting the nucleon-nucleus optical potential with low-momentum interactions is also inviting, given 
the improved status of Hartree-Fock as a starting point. 

Calculations with evolved cutoff-dependent 3N forces open the door to study 4N interactions. Pre- 
cision spectroscopy and low-energy reaction theory will require some control and assessment of the 
size of 4N forces. We note that the density dependence of 4N forces is not known for any low-energy 
nuclear Hamiltonian. Various approximations for evolving 4N interactions are suggested by analyzing 
the diagrammatic form of the SRG |137] . 

While the SRG is currently formulated in terms of coupled differential equations, experience with 
RG evolution for Viowfc suggest that an integral form may be more robust numerically. A possible 
starting point is the formalism developed in Ref. |257j . There are also many open questions concerning 
the relationship of the RG techniques discussed here to the wide range of RG technology applied to 
condensed matter and high-energy problems. For example, what are transformations that convert 
energy to momentum dependencies? And is there a path integral formulation of the SRG? 

The low-momentum RG technology that has been developed can also be usefully applied directly 
to EFT. Even though the widely available chiral EFT potentials have cutoffs as low as 500 MeV, these 
are not sharp cutoffs and significant high momentum strength is still present. As shown in Ref. [6], an 
exponential regulator with riexp ~ 8 leads to weaker distortions than the conventional value of riexp = 3 
used in the N^LO NN potentials [201 Hi]. At still lower cutoffs, an SRG- type regulator may have 

Calculations with low-momentum NN interactions are useful to establish convergence properties and to use the cutoff 
dependence to learn about the omitted forces, or in certain applications where 3N contributions are small, such as 
low-density neutron matter. 
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additional advantages. Experience with lower momentum regulators could be merged with expertise 
on fitting EFT low-energy constants to optimize the EFT at lower (softer) cutoffs for testing power 
counting as well as for applications to few- and many-body methods. This will help resolve questions 
about the utility of Viowk interactions compared to unevolved chiral EFT with lower cutoffs and possibly 
facilitate the use of error plots for the EFT [56] . 




For infinite matter, a power counting is needed for low-momentum interactions that specifies con- 
trolled many-body truncations at finite density. The Bethe-Brueckner-Goldstone expansion, which 
involves a summation of particle-particle ladders into the G matrix, shows (for potentials with strong 
short-range repulsion) that a perturbative expansion in terms of the G matrix does not work and so 
further summations are required. (The summation in this case is organized by the hole-line expan- 
sion.) Indications are that both summations are unnecessary with low-momentum interactions, but a 
quantitative power counting (for example, using the momentum dependence of the interactions together 
with phase-space constraints) needs to be demonstrated. This will involve understanding the evolution 
of particle-hole physics and what are the non-pertubative aspects of low-momentum interactions for 
nuclear structure. 

Finally, there are also opportunities for applications beyond ordinary nuclei. These include low- 
momentum interactions between baryons that carry strangeness for descriptions of hypernuclei (see, 
for example, Refs. |12 1(12581 - 1260] ). Renormalization group methods are already applied widely in con- 
densed matter physics, but the special experience from nuclei may be of advantage. The apphcation to 
microscopic inter-atomic potentials is an example that may prove fruitful. 

In closing, we reiterate that methods to derive low-momentum interactions not based (at least 
explicitly) on RG evolution are being pursued with promising results and prospects. The future interplay 
of the different approaches will be fruitful. When low-momentum interactions are coupled with the 
advances of many-body methods, the steady improvements in computational algorithms and power, 
and the on-going developments in chiral EFT and lattice QCD, the dream of microscopic calculations 
connecting QCD to the full range of nuclei and their reactions is within reach. 
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